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Abstract 

Nuclear weak interaction rates for /p-shell nuclei in stellar matter and the associated energy losses 
are calculated using a modified form of proton-neutron quasiparticle RPA model with separable Gamow- 
Teller forces. The stellar weak rates are calculated over a wide range of densities (10< pY e (gcm~ 3 ) < 
10 11 ) and temperatures (10 7 < T(K) < 30 x 10 9 ). This is the first ever extensive compilation of weak 
interaction rates in stellar matter calculated over a wide temperature-density grid and over a larger 
mass range. The calculated capture and decay rates take into consideration the latest experimental 
energy levels and ft value compilations. We take into consideration for the first time the effect of 
particle emission processes from excited states in the calculations of stellar rates which considerably 
affect our calculated rates. We have calculated stellar weak interaction rates for a total of 619 nuclei 
in the mass range A = 40 to 100. These also include many important neutron-rich nuclei which play 
an important role in the evolution process of stellar collapse. In a series of papers, starting from this 
one, we will be presenting our results on an abbreviated scale of temperature and density. This paper 
contains the stellar weak rates in the mass range 40 to 60. 
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1 INTRODUCTION 

The weak interaction has several crucial effects in 
the course of development of a star. It initiates the 
gravitational collapse of the core of a massive star 
triggering a supernova explosion, plays a key role in 
neutronisation of the core material via electron cap- 
ture by free protons and by nuclei and effects the 
formation of heavy elements above iron via the r- 
process at the final stage of the supernova explosion 
(including the so-called cosmochronometers which 
provide information about the age of the Galaxy 
and of the universe). The weak interaction also 
largely determines the mass of the core, and thus 
the strength and fate of the shock wave formed by 
the supernova explosion. 

The (3 decay properties of neutron-rich nuclei are 
a prerequisite for a better understanding of the r- 
process. The r-process theory was put forward by 
U more than four decades ago. Both the el- 
ement distribution on the r-path, and the result- 
ing final distribution of stable elements are highly 
sensitive to the decay properties of the neutron- 
rich nuclei involved in the process. This was also 
stressed by Klapdor || (see also and references 
therein). There are about 6000 nuclei between the 
/3 stability line and the neutron drip line. Most of 
these nuclei cannot be produced in terrestrial lab- 
oratories and one has to rely on theoretical extrap- 
olations in respect of beta decay properties. The 
microscopic calculations of weak interaction rates 
H ||, 0, D led to a better understanding of the r- 
process. However there was a need to go to domains 
of high temperature and density scales where the 
weak interaction rates are of decisive importance 
in studies of the stellar and galactic evolution pro- 
cesses and nucleosynthesis calculations. 

During the collapse to a neutron star or to a black 
hole, observation of the neutrino emission provides 
the only possibility for direct diagnosis of the events 
in the interior of the star. The neutrino spectrum 
with its flavors, its luminosities and its time depen- 
dence reflects the entire internal evolution and thus 
provides a simultaneous test for theories of stel- 
lar collapse and of the formation of neutron stars 
and black holes. The formation of neutron stars is 
one of the few instances where the weak interaction 
plays a crucial role in a macroscopic event in na- 
ture. It is generally assumed that massive stars in 
the region M > 8M Q (M Q denotes the solar mass) 



end in supernova explosions (of Type II). 

If the mass of the Fe core exceeds the so-called 
Chandrasekhar limit || 

M ch = 1.45(2T e ) 2 M (1) 

(Y e is the number of electrons present per nucleon) 
the pressure of the relativistic degenerate electron 
gas can no longer resist the force of gravity, and the 
core becomes unstable. The reason for the collapse 
is the photodisintegration of iron group nuclei and 
the onset of electron capture by free protons and 
nuclei |l(| which is made possible by the increase 
in the Fermi energy of the degenerate electron gas. 
The electron capture rates determine the initial dy- 
namics of the collapse and also, via Eq. (1), the size 
of the collapsing core-and therewith, the fate of the 
shock wave released later. These weak rates, ener- 
getically forbidden under earth conditions, increase 
by many orders of magnitude, due to the thermal 
population of parent excited states which also con- 
tribute to the overall rates. 

Prior to this work, the most extensive calcula- 
tions of stellar weak rates over a wide range of tem- 
perature (10 7 < T(K) < 10 11 ) and density (10< 
pY e (gem -3 ) < 10 11 ) was done by Fuller, Fowler 
and Newman (hereafter refered to as FFN). 
They calculated stellar electron and positron emis- 
sion rates and continuum electron and positron 
capture rates, as well as the associated neutrino 
energy loss rates for 226 nuclei with masses be- 
tween A — 21 and 60. Measured nuclear level infor- 
mation and matrix elements available at that time 
were used and unmeasured matrix elements for al- 
lowed transitions were assigned an average value 
of logft =5. To complete the FFN rate estimate, 
the Gamow- Teller contribution to the rate was pa- 
rameterized on the basis of the independent par- 
ticle model and supplemented by a contribution 
simulating low-lying transitions. The FFN rates 
were then updated and extended to heavier nuclei 
Q . These authors also considered the quenching 
of the Gamow- Teller strength by reducing the in- 
dependent particle estimate for the Gamow- Teller 
resonance contribution by a common factor of two. 
However there was a need for an improved theoreti- 
cal description of the problem [|l3| . The calculation 
of the weak interaction mediated reactions under 
the astrophysically reliant conditions requires the 
knowledge of not only the total strength in a given 
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direction (/3 _ -decay or electron capture) but also 
that of the strength distributed in nuclear excita- 
tion energy. Attempts have been made to obtain 
these distributions theoretically from shell model 
calculations. However, the number of shell model 
basis states can get very large for mid-/p nuclei, 
even at low energies. Hence the direct method 
of obtaining Gamow- Teller Giant Resonance dis- 
tribution from shell model calculations using full 
OTilu basis is computationally a formidable task even 
for a few nuclei. Several attempts to calculate the 
Gamow- Teller distribution from direct shell model 
calculations using a truncated shell model basis 
space have also been made (eg. Jl3)) which because 
of the above reasons, despite substantial computa- 
tional efforts, are as yet approximate. Shell model 
Monte Carlo technique [ll| was then used to 
further refine the calculations in the /p-shell. Re- 
cently KB3 interaction |l6j was used to perform cal- 
culations of electron capture and (3 decay rates for a 
few nuclei by shell model diagonalization in the fp- 
shcll 0^,^. These authors calculated the su- 
pernova electron capture and beta-decay rates for a 
few nuclei. In any case the calculation of the weak 
interaction mediated reaction rates for a moder- 
ately large number of nuclei required in astrophys- 
ical situations requires a straightforward and com- 
putationally manageable approach. 

Nuclei with A > 60 were in the beginning not in- 
cluded in the stellar evolution calculations because 
nuclei with A > 60 were not thought to become 
abundant enough to be important. It was shown 
later that these large nuclei also have important 
effects on the physics of core collapse [poj . These 
authors emphasized that calculations of presuper- 
nova evolution generate cores that are so neutron- 
rich (Y e reaches 0.42 and lower) that nuclei more 
massive than A = 60 must be considered, thus re- 
quiring an extension of the FFN rate tables. Soon 
work was started on calculation of (3 decay for nu- 
clei with A > 60 |2l| . These authors reported calcu- 
lation of (3 decay for 1 1 nuclei using an average beta 
strength function for typical presupernova matter 
density (p = 3 x 10 7 to 3 x 10 9 g cm ) and tem- 
perature (T = (2 to 5) x 10 9 K). Later a search 
was made for important weak interaction nuclei in 
presupernova evolution [|l2| . One of the difficulties 
faced by these authors in performing their study 
was the lack of enough (3 decay and electron cap- 
ture rates for nuclei with A > 60. The authors then 



had to made an estimate of these rates and approx- 
imate rates were derived. Since then there was a 
need to have microscopic stellar rates for /p-shell 
nuclei. 

This work is the first ever extensive calculation 
of stellar weak rates in the fp- and fpg-sheW nu- 
clei ranging from A = 40 to 100. A total of 619 
nuclei were considered for the calculation of stel- 
lar weak rates. These include also proton-rich and 
neutron-rich nuclei. The (3 decay of neutron-rich 
nuclei prevents the decrease of the overall electron 
fraction in the stellar core Nuclei with A > 

60 may make significant contributions to the elec- 
tron capture and (3 decay rates in the last stages 
of stellar evolution [^0] . In this work twelve differ- 
ent weak rates have been calculated for each parent 
nucleus. These include e ± -capture rates, /3 ± -decay 
rates, neutrino (anti-neutrino) energy loss rates, 
energies of beta-delayed proton (neutron) and the 
probabilities of these beta-delayed particle emission 
processes. The good comparison of pn-QRPA to 
experimental data |§| encouraged us to use pn- 
QRPA theory to calculate stellar weak rates. Fur- 
ther the pn-QRPA theory is supposed to work even 
better for /p-shell nuclei where the number of nu- 
cleon is fairly large. 

Section 2 presents the formalism of our rate cal- 
culations for various nuclear processes. The refer- 
ences used for accumulating experimental data and 
their incorporation to our calculation is treated in 
Section 3. Results and discussions are introduced 
in Section 4. Here we also compare our results to 
earlier compilations. We discuss the reliability of 
the calculated rates using the pn-QRPA model in 
Section 5. Section 6 finally summarizes our results. 

2 FORMALISM 

2.1 Assumptions 

The following main assumptions are made: 

(1) Only allowed Gamow- Teller and superallowcd 
Fermi transitions are calculated. It is assumed that 
contributions from forbidden transitions are rela- 
tively negligible. 

(2) It is assumed that the temperature is high 
enough to ionize the atoms completely. The elec- 
trons are not bound anymore to the nucleus and 
obey the Fermi-Dirac distribution (see also the dis- 
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cussions in pj|). At high temperatures (kT > 
1 MeV) , positrons appear via electron-positron pair 
creation, and the positrons follow the same energy 
distribution function as the electrons. 

(3) Neutrinos and antineutrinos escape freely 
from the interior of the star. Therefore, there are 
no (anti) neutrinos which block the emission of these 
particles in the capture or decay processes. Also, 
(anti)neutrino capture is not taken into account. 

(4) The distortion of electron (positron) wave- 
function due to the Coulomb interaction with a nu- 
cleus is represented by the Fermi function in phase 
space integrals. 

(5) Particle emissions from excited states are 
taken into account. This was neglected in the pre- 
vious calculations of weak interaction rates and is 
found to considerably affect the capture and de- 
cay rates. The effects of separation energy of pro- 
tons (S p ) and neutrons (S n ) on the stellar rates are 
considered. All excited states, with energy greater 
than S p or S n , are assumed to decay to the ground 
state by emission of protons or neutrons, respec- 
tively. Due to uncertainties in the calculation of 
energy levels and the effect of the Coulomb bar- 
rier for the case of proton emission, it is assumed 
that particles are emitted at excited energies 1 MeV 
higher (called E cr u throughout this section) than 
the particle decay channel. 

(6) For the parent nuclei, the cut-off excitation 
energy is set at E cr u . All excited states in daughter 
nuclei, with excitation energy (Ej < E cr n) decay 
directly to the ground state through 7 transitions. 
All excited energy states (Ej > E cr i t ) decay by 
emitting particles (protons or neutrons) with in- 
creasing kinetic energies. It is further assumed that 
either protons (if S p < S n ) or neutrons (if S n < S p ) 
are emitted from these excited states. The (3 decay 
of a possible isomer is not taken into account. 



where (ft)ij is related to the reduced transition 
probability By of the nuclear transition by 



(3) 



The D appearing in Eq. (3) is a compound expres- 
sion of physical constants, 



D = 



2ln2h 7 n 3 

2 K A~ ' 



and, 



Bij = B(F) ij + (g A /g v ) 2 B(GT) l 



(4) 



(5) 



where B(F) and B(GT) are reduced transition 
probabilities of the Fermi and Gamow- Teller tran- 
sitions respectively, 



B(F) ij 



2j^K'llE4ll<>l a , < 6 > 

k 

= ^- I |< i || II *>r, (7) 

k 

In Eq. (7), a k is the spin operator and t± stands 
for the isospin raising and lowering operator. The 
value of D=6295 s is adopted and the ratio of the 
axial- vector (g A ) to the vector (<jy) coupling con- 
stant is taken as 1.254. 

2.2.1 Phase space integrals 

The phase space integral (/y) is an integral over 
total energy, 



fii = J ™Vw 2 - l(w m -w) 2 F(±Z,w)(l-G T )dw, 

(8) 

for electron (upper signs) or positron (lower signs) 
emission, or 



2.2 Weak Decay in Stellar Matter 

The weak decay rate from the ith state of the par- 
ent to the jth state of the daughter nucleus is given 

Ji 3 {T lPl E f ) 



Xij — ln2- 



(/*)< 



(2) 



throughout subsection 2.2 natural units (h = c = m e = 
1) are adopted, unless otherwise stated, where m e is the 
electron mass. 



fa 



w\Jw 2 - l(w m + w) 2 F(±Z, w)G T dw, 



(9) 

for continuum positron (lower signs) or electron 
(upper signs) capture. 

In Eqs. (8) and (9), w is the total kinetic en- 
ergy of the electron including its rest mass, wi is 
the total capture threshold energy (rest+kinetic) 
for positron (or electron) capture. One should 
note that if the corresponding electron (or positron) 
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emission total energy, w m , is greater than -1, then 
wi = 1, and if it is less than or equal to 1, then 
wi =| w m |. w m is the total /3-decay energy, 



m d + E l - E j 



(10) 



where m p and Ei are mass and excitation energies 
of the parent nucleus, and and Ej of the daugh- 
ter nucleus, respectively. 

G+ and G_ are the positron and electron distri- 
bution functions, respectively. Assuming that the 
electrons are not in a bound state, these are the 
Fermi-Dirac distribution functions, 



G- 



[exp{ 



E-E 



kT 



) + !]■ 



G+ = [exp( 



E + 2 + E f . 
kT 



(11) 



(12) 



Here E = (w — 1) is the kinetic energy of the elec- 
trons, Ef is the Fermi energy of the electrons, T is 
the temperature, and k is the Boltzmann constant. 

In the calculations, the inhibition of the final neu- 
trino phase space is never appreciable enough that 
neutrino (or anti-neutrino) distribution functions 
had to be taken into consideration. F(±Z, w) are 
the Fermi functions and are calculated according to 
the procedure adopted by Gove and Martin Q . 

The number density of electrons associated with 
protons and nuclei is pY e NA, where p is the baryon 
density, Y e is the ratio of electron number to the 
baryon number, and Na is the Avogadro number. 



pYe 



1 



TT 2 N A h 



(G_ - G+)p 2 d P , (13) 



where p — (w 2 — l) 1 / 2 is the electron or positron 
momentum, and Eq. (13) has the units of moles 
cm . This equation is used for an iterative calcu- 
lation of Fermi energies for selected values of pY e 
and T. 

There is a finite probability of occupation of par- 
ent excited states in the stellar environment as a re- 
sult of the high temperature in the interior of mas- 
sive stars. Weak interaction rates then also have a 
finite contribution from these excited states. The 
occupation probability of a state i is calculated on 
the assumption of thermal equilibrium, 



P 



(2J, + l)exp(-Ei/kT) 
J2 i=1 (2Ji + l)exp(-Ei/kTy 



(14) 



where Ji and Ei are the angular momentum and 
excitation energy of the state i, respectively. 

Unfortunately one cannot calculate the J^'s in 
QRPA theory and hence Eq. (14) is modified as 



P 



exp(—E{/kT) 



(15) 



This approximation is a compromise and can be 
justified when one takes into consideration the un- 
certainty in the calculation of Ei which easily over- 
sheds the uncertainty in calculating the values of 
Ji in the above Eq. (14). 

The rate per unit time per nucleus for any weak 
process is then given by 



A 



(16) 



The summation over all initial and final states is 
carried out until satisfactory convergence in the 
rate calculations is achieved. 

2.2.2 Other calculated rates 

The neutrino energy loss rates are calculated using 
the same formalism except that the phase space 
integral is replaced by 



fij = J Ww 2 - l(w m -w) 3 F(±Z,w)(l-G T )dw, 

(17) 

for electron (upper signs) or positron (lower signs) 
emission, or by 



f" 



waJ w 2 — l(io TO + w) 3 F(±, w)G^:dw : 



(18) 

for continuum positron (lower signs) or electron 
(upper signs) capture. 

The Xij in Eq. (16) in this case is the sum of 
the electron capture and positron decay rates, or 
the sum of the positron capture and electron decay 
rates, for the transition i — > j. 

The proton energy rate from the daughter nu- 
cleus is calculated, whenever S p < S n , by 



A P — ^ Pi\ij(Ej — E cr it), 



(19) 



for all Ej > E cr it, whereas for all Ej < E cr u one 
calculates the 7 heating rate, 



A 7 



2_^Pi\ijEj. 



(20) 
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If on the other hand, S n < S p , then the neutron 
energy rate from the daughter nucleus is calculated 
using 

\ n = Y^Pi\j{E j - E crit ), (21) 

ij 

for all Ej > E crit , and for all Ej < E crit the 7 
heating rate is calculated as in Eq. (20). 

The probability of /3-delayed proton (neutron) 
emission is calculated by 

pp{n) = ( 22 ) 

where jf are states in the daughter nucleus for 
which Ej, > Ecru. In Eqs. [(19)-(22)] is 
the sum of the electron capture and positron decay 
rates, or the sum of positron capture and electron 
decay rates, for the transition i — > j(jf). 

2.3 Gamow-Teller transitions 

The parent excited states can be constructed as 
phonon-correlated multi-quasi-particle states. An 
extension of the pn-QRPA model is straight for- 
ward to the Gamow-Teller transitions from nuclear 
excited states. Low-lying excited states of a nucleus 
are obtained by one-proton (or one-neutron) exci- 
tations. They are described, in the quasi-particle 
(q.p.) picture, by adding two-proton (two-neutron) 
q.p.'s to the ground state j24|. The transition am- 
plitudes between the multi-quasi-particle states can 
be reduced to those of single-quasi-particle states. 
Details of the formalism can be found in j25| . 

In the calculation of the Gamow-Teller transi- 
tions a quenching of the transitions was not ex- 
plicitly taken into account. The quenching of the 
Gamow-Teller strength cannot be a constant renor- 
malization of the axial vector current [^6[ 0], We 
did not choose a global quenching factor (as done 
in many shell model calculations) for the following 
reasons. In order to reproduce the Gamow-Teller 
strength the Gamow-Teller interaction strength pa- 
rameters x (f° r particle-hole interactions) and k 
(for particle-particle interactions) were adjusted in 
our calculations. These values were deduced from a 
fit to experimental half-lives and for every isotopic 
chain fixed values of \ an d k allowed to deduce 
a locally best value of x an d K, respectively. With 
the large model space (up to 7 major shells) consid- 
ered in this QRPA calculation and the fine tuning 



of the Gamow-Teller strength parameter, an extra 
quenching factor on one hand might slightly im- 
prove the experimental half-lives for some (3 + de- 
cays of very proton-rich nuclei but the overall com- 
parison with measured rates of nuclei would remain 
more or less unaltered (see [§2|). 

2.4 Fermi transitions 

The Fermi operator is independent of space and 
spin [Eq. (6)], and as a result, the Fermi strength 
is concentrated in a very narrow resonance cen- 
tered around the isobaric analog state (IAS) for the 
ground and excited states (The IAS is generated by 
operating on the associated parent states with the 
isospin raising or lowering operator: 

T±=X)i±(i), 

i 

where the sum is over the nucleons. This operator 
commutes with all parts of the nuclear Hamiltonian 
except for the Coulomb part). The superallowed 
Fermi transitions were assumed to be concentrated 
in the IAS of the parent state. The Fermi matrix 
element depends only on the nuclear isospin, T, 
and it's projection T z = (Z — N)/2 for the parent 
and daughter nucleus. The energy of the IAS is 
calculated according to the prescription given in 
pages 111-112] whereas the reduced transition 
probability is given by 

B(F) = T(T+l)-T zi T zf , (23) 

where T Z i and T z f are the third components of the 
isospin of initial and final analogue states. 

3 INCORPORATION OF 
EXPERIMENTAL DATA 

Table (1) shows the references from which the latest 
experimental excitation energies and log ft values 
are taken (NP — > Nuclear Physics and NDS — > Nu- 
clear Data Sheets). For those nuclei where more 
than one references were given, the latest data was 
taken into account. 

The calculated excitation energies were replaced 
with the measured ones when they were within 0.5 
MeV of each other. The log ft value of this energy 
level was also then replaced by the measured one. 
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Very low lying states were inserted in the calcula- 
tions together with their log ft values if the theory 
was missing them. Inverse and mirror transitions 
were also taken into consideration. If there ap- 
peared a level in experimental compilations without 
definite spin and parity assignment, no theoretical 
levels were replaced with experimental levels be- 
yond this excitation energy nor were they inserted 
in the calculations. 

The Q-value of each transition as well as S p and 
S n of each nucleus were derived using the experi- 
mental mass compilation of Audi et al. p8| . For 
nuclei where Q failed to give mass defects Moller 
and Nix ]2^] and Myers and Swiatecki were 
used to derive the corresponding energies necessary 
for the calculations. The two different theoretical 
mass formulae used reflect the sensitivity of stellar 
rates on input masses. Work is in progress to com- 
pare the results of using these different mass models 
on the rate calculations for the case of neutron-rich 
and proton-rich nuclei. In (3^] the stellar rates of 
neutron-rich nuclei using the mass model of Moller 
and Nix (29| are presented. Stellar rates using mass 
formula of Myers and Swiatecki |50| , for the case of 
neutron-rich and proton-rich nuclei, are given in 
this paper on an abbreviated scale of density. The 
complete rate table on magnetic tapes can be re- 
quested from the authors. 



4 RESULTS AND DISCUS- 
SION 

For nuclei in the mass region A = 40 to 60, some 
work has been done previously regarding the cal- 
culation of stellar weak rates. The calculations in- 
clude the work of FFN (which to-date is used in 
many simulations for example the KEPLER stel- 
lar evolution code ]3S|) and the recently reported 
shell model diagonalization approach calculations 
jl7], [l8|, [l9|]. The shell model diagonalization ap- 
proach is used to calculate the stellar electron cap- 
ture and beta decay rates of a few nuclei. Our cal- 
culations are compared with the FFN and the shell 
model calculations (reported at the time of writing 
this paper) 0, [jj, . 



4.1 Comparison with the FFN calcu- 
lation 

Ten different cases of electron capture and beta- 
decay nuclei were considered for comparison with 
earlier calculations of FFN. The idea was to chose 
the most important weak interaction nuclei which 
play a key role in the evolution process of the cores. 
For these the nuclei in Table 25 and Table 26 of 
|l2| were considered. In these tables the authors 
have presented a list of 90 electron capture nuclei 
and 71 beta-decay nuclei averaged throughout the 
stellar trajectory for 0.40 < Y e < 0.5. These were 
the nuclei which, according to the calculations of 
the authors, affected Y e (rate of change of Y e ) the 
most in the presupernova evolution. These lists 
also contain nuclei with A > 60. The top 10 nuclei 
with masses 40 to 60 in both lists were chosen for 
comparison with this work. 

The QRPA rates are generally suppressed as 
compared to the rates of FFN. The QRPA rates 
are suppressed up to more than five orders of mag- 
nitude. There are two main reasons for this sup- 
pression. FFN does not take into effect the pro- 
cess of particle emission from excited states. The 
method for truncation of a state-density integral 
for calculation of the nuclear partition function was 
criticized by Tubbs and Koonin p3). The single 
particle level density, g(e), which is used for the 
calculation of the partition function, was discussed 
in detail in In (3^] it is argued that the no- 

truncation method, prescribed by J33|, is associated 
with the presence of very high angular momentum 
states. However in the Gamow- Teller transitions, 
considered in our work, only low angular momen- 
tum states are considered. Thielemann and collab- 
orators |35) also accept the consideration of ener- 
gies close to the particle separation thresholds for 
the calculation of level densities for astrophysical 
purposes. A choice of particle threshold decay as 
the cutoff parameter for parent excitation energy 
seems to be a reasonable choice as also discussed 
earlier by Fowler, Engelbrecht and Woosley p6| . 
This work takes into consideration for the first time 
the particle emission processes, which constrain the 
parent excitation energies. At any Ei, higher than 
the minimum of S p and S n (after accounting for 
the effective Coulomb barrier which prevents a pro- 
ton from being promptly emitted, and, the uncer- 
tainty in calculation of energy levels), the nucleus 
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Figure 1: Comparison of the QRPA electron capture (EC) 
rates (this work) with those of FFN. Solid lines represent 
the QRPA electron capture rates calculated in this work, 
while broken lines represent the electron capture rates of 
FFN. Log(T) is the log of temperature in units of Kelvin 
and log(EC) represents the log of electron capture rates in 
units of sec -1 . The upper graph is plotted at pY t 
cm -3 and the lower graph is plotted at pY e 



10 3 g 
10 11 g cm 3 . 



will emit protons or neutrons (as the case may be) 
instead of continuing to undergo capture process. 
The Ei considered in this work, as a result, are 
considerably lower as compared to those of FFN. 
Secondly for odd-A nuclei FFN places the centroid 
of the Gamow- Teller strength at too low excita- 
tion energies (discussed also in Jll], ) . Their 
rates are thus somewhat overestimated. The top 
10 electron capture nuclei are, in order of decreas- 
ing importance, 60 Co, 55 Co, 56 Ni, 57 Co, 55 Fe, 59 Co, 
54 Mn, 53 Mn, 54 Fe and 56 Fe in the mass range 40 to 
60 @. 

Fig. 1 shows the comparison of the QRPA elec- 
tron capture rates with those of FFN for the most 
important nucleus 60 Co. The QRPA rate is sup- 
pressed by five orders of magnitude at log(T) = 
9 and pY e = 10 3 g cm~ 3 . The suppression fac- 
tor decreases as density increases to pY e — 10 11 g 
cm~ 3 . FFN adopts the so-called Brink's hypothe- 
sis in their calculations. This hypothesis assumes 
that the Gamow- Teller strength distribution on ex- 
cited states is the same as for the ground state, only 
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Figure 2: Same as Fig. 1 but for electron capture by 56 Fc. 



shifted by the excitation energy of the state. This 
hypothesis is used because no experimental data is 
available for the Gamow- Teller strength distribu- 
tions from excited states and FFN did not employ 
any microscopic theory to calculate the Gamow- 
Teller strength functions from excited states. The 
pn-QRPA theory calculates stronger Gamow- Teller 
strength distribution from these excited states com- 
pared to those assumed using Brink's hypothesis 
and hence the suppression is reduced to a factor 
three at log(T) = 9 and pY e = 10 11 g cm" 3 . The 
FFN rates increase faster at higher temperatures as 
their parent excitation energies are not constrained, 
as mentioned above. 

Fig. 2, Fig. 3 and Fig. 4 show the comparison of 
electron capture rates for 56 Fe, 54 Mn and 54 Fe, re- 
spectively. These figures follow similar trends and 
one sees that the suppression factor at log(T) = 9 



and pY e = 10 g cm~ 3 keeps on decreasing from 
2 orders of magnitude ( 56 Fe) to about a factor of 
three ( 54 Fe) . The enhancement in the QRPA rates 
continues and for 57 Co and 56 Ni the agreement 
seems to be very good at log(T) = 9 and pY e = 10 3 
g cm -3 (Fig. 5 and Fig. 6, respectively). The trends 
are similar at pY e — 10 11 g cm -3 . This trend con- 
tinues and one sees that for 59 Co the QRPA elec- 
tron capture rate is enhanced by a factor of 5 in 
comparison to FFN at the lower density (Fig. 7). 
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Figure 4: Same as Fig. 1 but for electron capture by 54 Fe. Figure 6: Same as Fig. 1 but for electron capture by 56 Ni. 
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Figure 7: 

59 Co. 



Same as Fig. 1 but for electron capture by 



Clearly here the pn-QRPA theory calculates strong 
Gamow- Teller distribution strength function from 
parent excited states. However one notes that at 
higher temperatures, log(T) > 9.5, the FFN rates 
are again enhanced in comparison to our rates. At 
high temperatures the probability of occupation of 
parent excited states (Ei) increases [Eq. (14)] and, 
since FFN considers higher parent excited states, 
their rates are enhanced at higher temperatures. 

Fig. 8, Fig. 9 and Fig. 10 show a comparison of 
electron capture rates for the nuclei 55 Fe, 53 Mn and 
55 Co, respectively. The QRPA rates keep on en- 
hancing at pY e = 10 3 g cm -3 and low temperatures. 
In Fig. 10 the QRPA rates are enhanced by up to an 
order of magnitude in comparison to the FFN rates. 
At high temperatures the effect of high Ei consid- 
ered by FFN dominates and hence FFN electron 
capture rates exceed the QRPA rates. Comparison 
of the electron capture rates at p Y e — 10 11 g cm' 3 
remains more or less the same throughout the con- 
sidered cases. 

The top 10 beta-decay nuclei are, according to 
JlH , in order of decreasing importance, 56 Co, 59 Fe, 
&Mn, 56 Mn, 57 Fe, 57 Mn, 60 Co, 54 Cr, 58 Co and 
60 Fe in the mass region 40 to 60. It is very in- 
teresting to note that the pn-QRPA theory reports 
vanishing decay rates for 56 Co. This is in sharp 
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Figure 8: Same as Fig. 1 but for electron capture by 55 Fe. 
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Figure 9: 

53 Mn. 



Same as Fig. 1 but for electron capture by 
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Figure 10: Same as Fig. 1 but for electron capture by 
55 Co. 
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Figure 11: Comparison of the QRPA beta-decay (BD) 
rates (this work) with those of FFN. Solid lines represent 
the QRPA beta-decay rates calculated in this work while 
broken lines represent the decay rates of FFN. Log(T) is the 
log of temperature in units of Kelvin and log(BD) represents 
the log of /3~ rates in units of sec -1 . The upper graph is 
plotted at pY e = 10 3 g cm -3 and the lower graph is plotted 
at pY e = 10 11 g cm" 3 . 



contrast to Table 26 of jl^] where 56 Co is ranked 
as the most important beta-decay nuclei causing 
the largest change in Y e . The Q-value for this de- 
cay is -2.136 MeV and S p is 5.849 MeV. One of 
the selection rules for the quasi-particle transition 
to occur in the pn-QRPA theory is that the trans- 
forming quasi-particles should belong to the same 
major oscillator shell, i.e the neutron which changes 
to a proton should belong to the same major os- 
cillator shell (and same for proton). Even though 
the decay channel opens below the particle decay 
channel no decay rates are obtained for 56 Co. In 
this case the above mentioned selection rule is sat- 
isfied only for high-lying major shells but then the 
parent excited state (which is the sum of the en- 
ergies of constituent quasi-particles) is above the 
particle decay energy. 56 Co can decay however to 
56 Ni via positron capture process at higher temper- 
atures. So the top 10 beta-decay nuclei taken for 
comparison with FFN are, in order of decreasing 
importance, 59 Fe, 54 Mn, 56 Mn, 57 Fe, 57 Mn, 60 Co, 
54 Cr, 58 Co, 60 Fe and 57 Cr. 

The beta decay rates calculated in this work are 
also, in general, suppressed in comparison to the 
calculations of FFN by up to 13 orders of magni- 
tude and more. 

Fig. 11 depicts the comparison of the QRPA 
beta-decay rates (this work) and those of FFN for 
the odd-odd cobalt isotope, 58 Co. One notes that 
the QRPA rates are very much suppressed in com- 
parison to FFN rates. In the figure the QRPA rate 
is suppressed by some 13 orders of magnitude at 
log(T) = 9 and pY e = 10 3 g cm" 3 . At log(T) = 9 
and pY e = 10 11 g cm" 3 , the QRPA rate is van- 
ishing. Beta-decay rates are sensitive functions of 
(Ei — Ej). Even for positive Q-values where mea- 
sured Gamow- Teller transitions might also exist, 
the phase space (= Q + Ei — Ej) at high temper- 
atures and densities might be zero if Ei is consid- 
erably less than Ej. In this work E^s are con- 
strained due to the particle emission process and 
no Ei above particle decay channel (after account- 
ing for the effective Coulomb barrier and the un- 
certainty in calculation of energy levels) are con- 
sidered. Especially at high densities, where the 
chemical potential of electrons increases outside the 
nuclei (thus impeding the decay process), this phe- 
nomenon is very obvious. By a mere increment of 
500 keV in Ei many of the QRPA rates at high den- 
sities can increase by some tens of orders of mag- 
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nitude. At pY e — 10 11 g cm" 3 , the phase space for 
calculation of beta-decay rates is negative resulting 
in vanishing decay rates till log(T) = 9. Only at 
higher temperatures (log(T) > 9) the phase space 
becomes positive and then the QRPA rates tend to 
FFN rates. 

Fig. 12, Fig. 13 and Fig. 14 depict the compar- 
ison of beta-decay rates for 54 Mn, 54 Cr and 60 Fe, 
respectively. At lower density of pY e = 10 3 g cm" 3 , 
one notes that at log(T) = 9 the suppression rate 
is reducing gradually. This trend continues and fi- 
nally in Fig. 19 (representing 57 Cr) the QRPA rate 
is suppressed by a mere factor of 3 at pY e = 10 3 
g cm" 3 , log(T) = 9. Fig. 20 (representing 57 Mn) 
then depicts a situation where the QRPA rates are 
enhanced in comparison with FFN rates at pY e = 
10 3 g cm" 3 till log(T) = 9.5. At higher tempera- 
tures FFN rates increase due to the phenomenon 
of considering higher E^s in their calculation. In 
all these figures one sees that the situation remains 
more or less the same at pY e = 10 11 g cm" 3 . 

The comparison of positron capture and positron 
decay rates is similar to the preceding comparison 
of electron capture and beta decay rates, respec- 
tively. In our table of rates, - 100 means that the 
rate is smaller than 10 -100 . It should be pointed 
out that the decay rates are sensitive functions 



Figure 13: Same as Fig. 11 but for beta decay by 54 Cr. 
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Figure 14: Same as Fig. 11 but for beta decay by 60 Fe. 
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Fi glire 16: Same as Fig. 11 but for beta decay by 56 Mn. Figure 18: Same as Fig. 11 but for beta decay by 60 Co. 
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Figure 19: Same as Fig. 11 but for beta decay by 57 Cr. 
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Figure 20: Same as Fig. 11 but for beta decay by 57 Mn. 



of the difference of parent and daughter energies 
(Ei — Ej) and just by a mere addition of around 
0.5 MeV, many of our -100 decay rates could rise 
to larger numbers. For all cases where experimen- 
tal data were available our results are in very good 
agreement with the results of FFN. We differ when 
no experimental data are present. In all such cases 
FFN uses either some parameterization or some ap- 
proximation for estimation of nuclear matrix ele- 
ments whereas the nuclear matrix elements are cal- 
culated in the pn-QRPA theory in our compilation. 



4.2 Comparison with the shell model 
diagonalization approach calcu- 
lation 

The shell model calculations describe well the nu- 
clear structure and, especially for the light nuclei 
(eg., ]3^]), agree very well with the experiment. 
The dimension of the many-body space, however, 
grows exponentially as the number of particle in- 
creases, making it difficult to apply shell model 
techniques to heavier nuclei. A comparison of half- 
lives calculated by the pn-QRPA model for light 
nuclei with shell model calculations can be seen in 
Figure 8 of p^] . Recently the shell model diagonal- 
ization approach (SMDA) has been used to calcu- 
late P decay and electron capture rates of some fp- 
shell nuclei |l§|l§. For the case of odd-odd nu- 
clei, no experimental information exists about the 
Gamow- Teller strength distribution. Comparison 
with the QRPA electron capture rates for odd-odd 
nuclei suggests that SMDA rates are suppressed. 
The density of states in odd-odd nuclei is quite 
high and doing calculations for a few Ei is unde- 
sirable. Further the shell model gives very weak 
transition strengths to the low-lying states in the 
daughter nucleus. Due to the model used, they 
also had to consider very few Ei in calculation of 
capture rates on odd-A and even-even nuclei (for 
56 Ni, they considered just the ground-state transi- 
tion). Consequently their rates are generally sup- 
pressed in comparison to the earlier calculations of 
FFN and |l2j . Similar conclusions can be drawn for 
the SMDA calculation of /3-decay rates. Here they 
considered Ei usually up to 1 MeV and "back res- 
onances" (the Gamow- Teller back resonances are 
states reached by the strong Gamow- Teller tran- 
sitions in the electron capture process built on the 
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Figure 21: Comparison of the Q 
(EC) rates (this work) with those of 
sent the QRPA electron capture rates oFthis work while bro 
ken lines represent the electron capture rates of Jl7j ]. Log(T) 
is the log of temperature in units of Kelvin and log(EC) rep- 
resents the log of electron capture rates in units of sec -1 . 
Aden represents the log of density in units of g cm -3 . 



Figure 22: 

55 Co. 



Same as Fig. 21 but for electron capture by 



§. These states 
) from daugh- 



ground and excited states, see [|38| 
were discussed earlier by Klapdor 
ter states of energy below 1 MeV. Back resonances 
from daughter states higher than 1 MeV were not 
considered. The calculations were done for the tem- 
perature range, T9 = 1-10 (T9 is temperature in 
units of 10 9 K) and density range, pr = 10-1000 
(p7 is the density in units of 10 7 g cm -3 ). Shell 
model calculations of electron capture rates were 
reported in [jl7| [l^] . A comparison of their electron 
capture rates with the QRPA electron capture rates 
is presented below (also see pof). 

The authors in |l7j calculated the Gamow- Teller 
strength distributions for the ground state and first 
excited states in 55 Co and 56 Ni. The authors quote 
a noticeably smaller rate than the usually adopted 
rate for 55 Co. The authors claim that in odd-A 
nuclei the Gamow- Teller strength resides at higher 
energies than assumed by FFN. As seen in Fig. 21 
and Fig. 22, the shell model rates are suppressed 
compared to the pn-QRPA rates by an average fac- 
tor of 65 at log(T) = 1. At log(T) = 10 the rates are 
suppressed on the average by a factor of 3. Due to 



the model limitations, the authors restricted them- 
selves to parent excited states of a few MeV. They, 
for example, went to an £!j of 2.565 MeV in the case 
of 55 Co whereas the QRPA rates considered Ei as 
high as 6.064 MeV (after which protons should be 
emitted from 55 Co). 

The story is somewhat different for 56 Ni. Fig. 23 
depicts the comparison. Till log(T) = 3 the shell 
model rates are suppressed. At log(T) = 1, eg., 
shell model rates are suppressed on the average by a 
factor of 3. The comparison at log(T) =3 is perfect 
and then shell model rates start increasing. As the 
authors argue, for the even-even case the centroid 
of the Gamow- Teller strength is placed at too high 
excitation energies in FFN and |l2|. Placement of 
the Gamow- Teller strength at a lower excitation 
energy can explain the sudden increase in their cal- 
culated rates in comparison with the QRPA rates 
with increasing temperature. 

Fig. 24 shows the comparison of the electron 
capture rates of the odd-odd cobalt isotopes, 
56,58,60 Co _ In thc nrst two graphs trie qrpa rates 

are enhanced because all possible excited states in 
parent nuclei were considered for calculation in this 
work which were below the particle decay chan- 
nel, whereas in the work of the first few ex- 
cited states were considered. In the last case at 
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Figure 23: Same as Fig. 21 but for electron capture by 
56 Ni. 
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Figure 24: Same as Fig. 21 but for electron capture by 
odd-odd cobalt isotopes. 



higher temperatures (log(T) > 3) our rates surpass 
the shell model rates. The probability of occupa- 
tion of excited states increases with increasing tem- 
perature [Eq. (15)]. Considerable enhancements to 
the rates come from transitions from these excited 
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Figure 25: Same as Fig. 21 but for electron capture by 
odd-odd manganese isotopes. 



Fig. 25 which shows a comparison of electron cap- 
ture rates for odd-odd manganese isotopes is simi- 
lar in explanation to Fig. 24. No experimental in- 
formation exists for the Gamow- Teller strength for 
odd-odd nuclei and one has to rely totally on the 
theoretical model, and the truncations might cause 
the suppression in their rates. Further for the case 
of odd-odd nuclei, the shell model gives very weak 
transition strength to the low lying states in the 
daughter nucleus fL8fl . 

The calculation of f3 decay rates using SMDA was 
discussed in |l]|. Here they tried to enhance their 
rates by using the back resonances from daughter 
states below 1 MeV energy. However no estimate of 
the energy was given to which the back resonances 
reach. Figures 26-28 depict a comparison of the (3 
decay rates. As can be seen (Fig. 26) the QRPA 
rates are suppressed by about one order of magni- 
tude. The same trend is seen in the comparison 
of odd-A nuclei (Fig. 27). The results support the 
statement in that experimental data and shell 
model studies tend to decrease the contribution of 
the back-resonance to the (3 decay rates for even- 
even and odd-A nuclei. One should note that at 
high densities and temperatures the QRPA rates 
surpass the SMDA rates. This shows the necessity 
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Figure 26: Comparison of the QRPA beta-decay (BD) 
rates (this work) with those of | JL9[ | . Solid lines represent the 
QRPA /3 rates of this work while broken lines represent the 
/3 rates of ]l9| , Log(T) is the log of temperature in units 
of Kelvin and log(BD) represents the log of /3 rates in units 
of sec -1 . Aden represents the log of density in units of g 
cm -3 . 



of considering transitions from high-lying states at 
high temperatures and densities. The QRPA ap- 
proach, used in the present calculation, allowed the 
use of a very large model space considering exci- 
tation of particles up to 7 major shells (SMDA, in 
comparison, considered 1 major shell). 

The difference in comparison is greatest for the 
case of odd-odd nuclei. Here the QRPA rates are 
suppressed by more than 7 orders of magnitude 
(Fig. 28). Our results do not support the conclu- 
sions drawn in JT9| that the (3 decay rate of odd-odd 
nuclei should be large compared to previous compi- 
lations. Table (2) and Table (3) show the compar- 
ison of the QRPA rates^SMDA rates, FFN rates 
and the calculation of [Q for electron capture and 
(3 decay rates, respectively. Table (2) and Table (3) 
are similar to Table I and Table II of [[l9| . 

The probabilities of /3-delayed proton (neutron) 
pp(«) emission at temperature 10 7 K and density 
pY e — 10 g cm -3 can be compared to terrestrial 
probabilities. We searched the literature for mea- 
sured probabilities to check the reliability of the 
present calculation. Unfortunately measured prob- 
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Figure 27: Same as Fig. 26 but for beta decay by odd-A 
nuclei. 
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Figure 28: Same as Fig. 26 but for beta decay by odd-odd 
nuclei. 



abilities are not always given; often it is only stated 
that particle emission processes are observed. We 
found that agreement with available experiment is 
satisfactory and that in all instances of reported 
particle emission our calculation also gives finite 
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probabilities pp("). (See [ET| for a discussion of 
comparison of such probabilities with experimental 
observations.) 

5 RELIABILITY OF THE 
pn-QRPA CALCULA- 
TIONS 

It should be pointed out that the uncertainty in 
calculations of stellar weak rates can be consider- 
ably large. Lack of experimental data deteriorates 
the situation. The "electron capture direction" can 
be explored experimentally by (n,p) experiments, 
whereas the "f3~ decay" side can be explored by 
the (p,n) reactions. There are a handful of other 
experiments which have been used by many theo- 
rists to shape the centroid and width of the Gamow- 
Teller strength. However these experimental data 
are not enough to completely explore the domain 
of nuclei which are interesting from astrophysical 
viewpoint. There exists, eg., no experimental in- 
formation about the Gamow- Teller strength distri- 
bution for odd-odd nuclei in the /p-shell. Further 
there is almost no experimental data concerning 
the Gamow- Teller strength distribution from par- 
ent excited states. In the stellar environment, at 
high temperatures and densities, there is a finite 
probability of occupation of parent excited states 
and transitions from these excited states are some- 
times many orders of magnitude higher than tran- 
sitions from the ground state. For most of the cal- 
culations then one uses the so-called Brink hypoth- 
esis which assumes that the Gamow- Teller strength 
distribution on excited states is the same as for the 
ground state, only shifted by the excitation energy 
of the state. This hypothesis is a very crude ap- 
proximation and sometimes can lead to misleading 
results. 

The weak interaction rates are calculated using 
Ay = -jjfijBij. (24) 

Here i represents the parent excites states and j the 
daughter's. The first factor is a constant and the 
second are phase space integrals which can be cal- 
culated relatively accurate. In this work the weak 
interaction rates are calculated using a microscopic 
theory. This theory (pn-QRPA) constructs parent 



and daughter excited states and also calculates the 
Gamow- Teller strength distribution among these 
states. In other words the Brink hypothesis is not 
employed in this calculation. However there are 
still uncertainties present. There are, eg., uncer- 
tainties present in the calculation of excited states. 
Mass defects are used as input parameters to cal- 
culate the Q-value of the reaction and S n and S p 
of parent and daughter nuclei. For certain neutron- 
rich and proton-rich nuclei, no experimental masses 
are present and one has to rely on theoretical mass 
models. The uncertainty in the Q-values must be 
viewed as the limiting factor for the prediction of 
weak interaction rates of unknown isotopes (see 
). Therefore this work incorporates the latest ex- 
perimental compilations of excitation energies and 
log ft values wherever present (see Table (1) and 
Section 3), to ensure the most reliable calculation 
of stellar weak rates allowed by the pn-QRPA the- 
ory. 

There, however, exists a check for the calcula- 
tions of these weak interaction rates. At the lowest 
temperature considered in this work (T9 = 0.01), 
excited parent states are not appreciably popu- 
lated, while at the lowest density considered in this 
work (pY e = 10 gem -3 ), the continuum electron 
density is quite low, and except for electron capture 
, the stellar rates should be close to the terrestrial 
values. The exception arises because only contin- 
uum electron capture and no bound states capture 
is calculated. There is a reasonable amount of data 
of measured half-lives available as compared to the 
measured Gamow- Teller strength functions. The 
idea was to check the calculation of the half-lives 
using the pn-QRPA theory with the measured ones 
and to also check the predictive power of the pn- 
QRPA theory (regarding the half-lives) with the 
recent measurements. We present such compar- 
isons for thousands of nuclei in this section to give 
a feeling how the pn-QRPA theory can perform in 
the domain of high temperature and density for the 
calculation of stellar rates of 719 nuclei considered 
in this project. We also present the comparison of 
our calculated total Gamow- Teller strength func- 
tion with those measured from (n,p) and (p, n) re- 
actions at the end of this section. 

A comparison between the terrestrial half-lives 
(calculated by ^]) and the stellar half-lives at 
temperature (T9 = 0.01) and density (pY e = 10 
gem -3 ) is given in Table (4) and Table (5). The 
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terrestrial half-lives are calculated using the mass 
formula of p{| whereas the stellar rates are pre- 
sented using the mass formulae of [^9| ||(|. The 
slight differences are attributed to the fact that in 
the present calculation a model space of 7 major 
shells was considered for all nuclei, while the cal- 
culation for terrestrial half-lives || ^ considered 
sometimes also smaller model spaces. 
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Figure 29: Predictive power of the pn-QRPA the- 
ory. Here the ratios of the calculated f3 + /(EC) 
decay half-lives to measured ones are shown as a 
function of mass numbers. The measured values 
were obtained after the calculations ||. 

The terrestrial half-lives are compared to the ex- 
perimental data in Table (6) and Table (7). In Ta- 
bles (6) and (7), N denotes the number of experi- 
mentally known half-lives shorter than the limit in 
the second line, n is the number (and percentage) 
of isotopes reproduced under the condition given 
in the first column, and x is the average deviation 
defined by 
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(25) 



i=i 



where 



rpcal /rpexp 
i l/2' 1 l/2 



_ rpexp /rpcal 
— J l/2 / J l/2 



if 



if 



rpcal »> rpexp 
1 l/2 — 1 l/2 



rpcal , rpexp 
J l/2 *- J l/2 ■ 



For example, the pn-QRPA reproduces 93% (75%) 
of all experimentally known half-lives shorter than 
1 minute for /3+/EC within a factor of 10 (2) with 
an average deviation of x = 1.72 (1.31) and 96% 
(82%) of all known /3 _ -decaying nuclei with half- 
lives less than a minute are reproduced within a 
factor of 10 (2) with an average deviation of x = 
1.67 (1.36). 
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QRPA-Calculation 
mass formula: Hilf et al. 



25 50 75 100 125 150 175 200 225 250 

atomic mass number A 

Figure 30: Predictive power of the pn-QRPA the- 
ory for proton-rich nuclei. Here the ratios of the 
calculated (3 + /(EC) decay half-lives (using mass 
formula of |76|) to measured ones are shown as a 
function of mass numbers. The measured values 
were obtained after the calculations p|. 



It can be seen from the tables that the model 
works better with increasing neutron excess (cor- 
responding to shorter half- lives), that is, with in- 
creasing distance from the line of stability. This is 
in agreement with the expectation, since forbidden 
transitions arc neglected in the calculation. This is 
also a promising feature with respect to the predic- 
tion of unknown half-lives, implying that the pre- 
dictions are made on the basis of a realistic physical 
model (see for example, Table B and Fig. 6 of ||). 

The predictive power of the pn-QRPA model is 
shown in Table (8) . Here a comparison of the mea- 
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QRPA-Calculation 
mass formula: Grooteetal, 
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QRPA-Calculation 
mass formula: HilfetaL 



25 SO 75 IOC 125 150 175 200 225 250 

atomic mass number A 



25 50 75 100 125 150 175 200 225 250 

atomic mass number A 



Figure 31: Same as Fig. 30 but here the mass for- 
mula of WTl was used for the calculated half- lives. 
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QRPA-Calculation 

mass foranik: Moeller and Nix 



25 50 75 100 125 150 175 200 225 250 

atomic mass number A 

Figure 32: Same as Fig. 30 but here the mass for- 
mula of |E9| was used for the calculated half-lives. 



Figure 33: Same as Fig. 30 but depicting the predic- 
tive power of the pn-QRPA theory for neutron-rich 
nuclei (f3~ decays). 
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QRPA-Calculation 
mass formula: Groote et al. 



25 50 75 100 125 



150 175 200 225 250 

atomic mass number A 



Figure 34: Same as Fig. 33 but here the mass for- 
mula of l|77|| was used for the calculated half- lives. 



sured f3 + /(EC) decay half-lives with those calcu- 
lated in the pn-QRPA theory is given. The Q- values were taken from [El 



All the calculated 



20 



is 10 



* 'it* 



• i 

< * * 



QRPA-Calculaiion 

mass formula: Moeller and Nix 



. i .... i i , 



25 50 75 100 125 150 175 IQQ 225 250 

atomic mass number A 

Figure 35: Same as Fig. 33 but here the mass for- 
mula of [E9fl was used for the calculated half- lives. 



rates shown were computed before their mea- 
surements given in the third column. Fig. 29 shows 
that most of the calculated half-lives are in very 
good agreement with the measured ones. 

Also the j3 + /(EC) decay half-lives of many 
proton-rich nuclei were measured after they were 
calculated in Since experimental masses were 
not available for these proton-rich nuclei at that 
time, Hirsch et al. calculated these half-lives us- 
ing three different mass formulae from [fT6f , |77| 
and (29). Table (9) shows the comparison of these 
calculated half-lives with the measured ones. Fig- 
ures 30-32 show the comparison with these later 
measured experimental data for the calculations of 
half-lives using the mass formulae of J76|, |77j and 
p9| , respectively. Again the comparison is very sat- 
isfactory. 

A search of new literature was also made for the 
measurement of neutron-rich nuclei which were not 
measured at the time of their calculations in 0] 
to check the predictive power of the pn-QRPA for 
neutron-rich nuclei. Table (10) shows the measured 
values and those calculated earlier by using mass 
formulae of f[% and Q whereas the compar- 
ison is depicted graphically in Figures 33-35. Once 
again the comparison is very satisfactory showing 



the good predictive powers of the pn-QRPA the- 
ory and putting weight into the reliability of our 
calculations. 

Finally we also present the comparison of our cal- 
culated Gamow- Teller strength function with those 
measured. As stated earlier, not much experimen- 
tal measurements of the Gamow- Teller strength 
function are present, (p, n) and (n,p) reactions can 
be used to measure the total B{GT)~ and B(GT) + 
strength function, respectively. One should be cau- 
tioned about the uncertainties present in these sorts 
of measurements and various energy cutoffs are 
used as a reasonable upper limit on the energy at 
which Gamow- Teller strength could be reliably re- 
lated to measured AL = cross sections. Table 11 
and Table 12 show the comparison of our calculated 
Gamow- Teller strength function with the measured 
ones for (n,p) and (p, n) reactions, respectively. In 
the comparisons we used the same energy cutoff 
as used in the respective experiments. We also in- 
clude in these tables the results of the large-scale 
shell model calculations of the total Gamow- Teller 
strength function |7^] where ever reported. One 
notes that our comparison with measured strength 
for 60 ' 62 ' 64 Ni isotopes is not in good agreement. 
Mass formula |^9| predicts these proton-magic iso- 
topes to be essentially spherical. However this 
is not the case when one measures the electric 
quadrupole transition probabilities (see discussions 
in S). This might can cause the enhancement in 
our calculated strength function (around a factor 
of 2). Low- lying transitions can contribute consid- 
erably in stellar rates due to phase space enhance- 
ments. We should like to point here again that in 
this compilation of stellar rates, experimental in- 
formations, derived from measured half-lives, were 
used, when ever available, and inserted if the pn- 
QRPA theory was missing them (see Section 3). 



6 SUMMARY 

In a series of paper beginning with a paper consid- 
ering the sd-shell nuclei we present extensive 
microscopic calculations of weak interaction rates 
in stellar environment and the first ever one for the 
/p-shell nuclei. The calculations are done for 619 
nuclei in the mass range A = 40 to 100. In this first 
paper of the /p-shell series we present the results 
for nuclei in the mass range A = 40 to 60. The 
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pn-QRPA theory which gives good results for the 
terrestrial rates also promises to be a good candi- 
date for calculations of stellar rates. 

We compared our results to the FFN calculations 
and to the recently reported shell model calcula- 
tions. We considered for the first time the effect 
of particle emission processes from parent excited 
states in these types of calculations. This leads to 
a general suppression of decay rates, and electron 
capture rates at high temperatures and densities, in 
comparison to the FFN rates. Our calculations are 
in very good agreement with the FFN calculations 
wherever experimental data are available. 

These calculations, including weak interaction 
rates of many important neutron-rich nuclei, are of 
prime importance for stellar and galactic evolution 
processes and nucleosynthesis calculations. 

In the next papers of this series we plan to 
present our results for A = 61 to 80 and A = 81 to 
100, respectively. 
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Table (1): References of experimental data incorporated in this calculation (NP — > Nuclear Physics, 
NDS -> Nuclear Data Sheets). 



Mass Number A 


Reference 


40 - 44 


NP A521,l (1990) 


45 


NDS 65,1 (1992), NDS 40,149 (1983) 


46 


NDS 68,271 (1993), NDS 49,237 (1986) 


A 7 

4 ( 


1YTT\C i A 1 /i nnt;^ attic i /'iriQfiA 
INDb /4,1 (lyyoj, INJJb 48,1 (lyooj 


48 


NDS 68,1 (1993), NDS 45,557 (1985) 


49 


NDS 76,191 (1995), NDS 48,569 (1986) 


50 


NDS 75,1 (1995) 


51 


NDS 81,183 (1997) 


52 


NDS 71,659 (1994), NDS 58,677 (1989) 


53 


NDS 61,47 (1990), NDS 43,481 (1984) 


54 


NDS 68,887 (1993), NDS 50,255 (1987) 


55 


NDS 64,723 (1991) 


56 


NDS 67,523 (1992), NDS 51,1 (1987) 


57 


NDS 67,195 (1992), NDS 47,1 (1986) 


58 


NDS 80,789 (1997) 


59 


NDS 69,733 (1993), NDS 39,641 (1983) 


60 


NDS 69,1 (1993), NDS 48,251 (1986) 
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Table (2): Comparison of the electron capture rates calculated by the pn-QRPA theory (this work) 
with previous works, pi is the density (in units of 10 7 g cm -3 ), Tg is the temperature (in units of 10 9 K). 
QRPA, SM, FFN and AFWH denote the rates calculated in this work, and Q, respectively. 
Exponents are given in parenthesis. All rates are given in units of s~ 4 . 



Nucleus 


f>7 


Tg 


QRPA 


SM 


FFN 


AFWH 


56 Ni 


4.32 


3.26 


9.9 


(-3) 


1.3 


(-2) 


7.4 


(-3) 


8.6 


(-3) 


54 Fe 


5.86 


3.40 


1.3 


(-5) 


4.2 


(-5) 


2.9 


(-4) 


3.1 


(-4) 




5.86 


3.40 


3.7 


(-4) 


8.1 


(-5) 


3.7 


(-4) 


6.3 


(-4) 


56p e 


10.7 


3.65 


1.1 


(-6) 


2.1 


(-6) 


1.0 


(-5) 


4.7 


(-7) 


bb Co 


4.32 


3.26 


8.0 


(-2) 


1.6 


(-3) 


8.4 


(-2) 


5.1 


(-2) 


57 Co 


5.86 


3.40 


1.6 


(-3) 


1.3 


(-4) 


1.9 


(-3) 


3.4 


(-3) 


55 Fe 


5.86 


3.40 


4.8 


(-3) 


1.9 


(-4) 


5.8 


(-3) 


3.8 


(-3) 




5.86 


3.40 


4.1 


(-3) 


4.7 


(-4) 


4.4 


(-3) 


4.4 


(-3) 


59 Co 


10.7 


3.65 


4.9 


(-4) 


7.8 


(-6) 


2.1 


(-4) 


2.1 


(-4) 


53 Mn 


10.7 


3.65 


1.4 


(-2) 


3.3 


(-4) 


3.8 


(-3) 


5.6 


(-3) 


5b Co 


5.86 


3.40 


3.3 


(-2) 


1.7 


(-3) 


6.9 


(-2) 


5.1 


(-2) 


54 Mn 


10.7 


3.65 


7.5 


(-4) 


3.1 


(-4) 


4.5 


(-3) 


1.1 


(-2) 


58 Co 


10.7 


3.65 


3.4 


(-3) 


3.5 


(-4) 


9.1 


(-3) 


2.1 


(-2) 


56 Mn 


33.0 


4.24 


1.1 


(-2) 


1.0 


(-4) 


4.1 


(-4) 


2.0 


(-3) 


60 Co 


33.0 


4.24 


2.0 


(-3) 


1.7 


(-4) 


1.1 


(-1) 


6.1 


(-2) 



Table (3): Comparison of the (3 decay rates calculated by the pn-QRPA theory (this work) with previous 
works. p7 is the density (in units of 10 7 g cm" 3 ), Tg is the temperature (in units of 10 9 K). QRPA, SM, 
FFN and AFWH denote the rates calculated in this work, pll, |llj and jj^ , respectively. Exponents 
are given in parenthesis. All rates are given in units of s . FFN did not calculate rates for nuclei with 
A > 60. 



Nucleus 


P7 


Tg 


QRPA 


SM 


FFN 


AFWH 


56p e 


5.86 


3.40 


5.5 (-13) 


3.9 (-11) 


2.3 (-10) 


5.9 (-11) 


54 Cr 


5.86 


3.40 


1.8 (-8) 


2.2 


(-7) 


2.2 (-5) 


1.5 


(-7) 


58 Fe 


10.7 


3.65 


7.3 (-9) 


5.2 


(-8) 


2.6 (-6) 


1.5 


(-7) 


eo Fe 


33.0 


4.24 


1.1 (-5) 


1.7 


(-4) 


4.6 (-3) 


1.0 


(-3) 


5 2 Ti 


33.0 


4.24 


1.8 (-5) 


1.3 


(-3) 


1.1 (-2) 


1.2 


(-4) 


59 Fe 


33.0 


4.24 


6.2 (-5) 


6.0 


(-5) 


6.3 (-3) 


5.3 


(-3) 


61 Fe 


33.0 


4.24 


4.2 (-3) 


1.7 


(-3) 




6.4 


(-2) 


61 Co 


33.0 


4.24 


2.1 (-5) 


1.6 


(-4) 




9.3 


(-4) 


63 Co 


33.0 


4.24 


3.8 (-2) 


1.6 


(-2) 




1.4 


(-2) 


59 Mn 


220 


5.39 


1.1 (-2) 


2.2 


(-2) 


7.2 (-1) 


1.4 


(-1) 


5 «Co 


4.32 


3.26 


2.7 (-11) 


2.7 


(-6) 


1.2 (-6) 


3.8 


(-6) 


54 Mn 


5.86 


3.40 


1.8 (-10) 


2.7 


(-6) 


1.6 (-6) 


7.5 


(-6) 


56 Mn 


10.7 


3.65 


5.7 (-6) 


3.4 


(-3) 


3.0 (-3) 


9.1 


(-3) 


60 Co 


10.7 


3.65 


8.3 (-7) 


6.6 


(-4) 


1.4 (-3) 


3.4 


(-3) 


50 Sc 


33.0 


4.24 


6.6 (-4) 


1.2 


(-2) 


2.8 (-2) 


1.8 


(-1) 
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Comparison of the terrestrial and stellar half-lives at lowest temperature and density 

Table (4): Comparison of calculated j3 + /{EC) decay half-lives at temperature (T9 = 0.01) and density 
(pY e — 10 gem -3 ), with the calculation of terrestrial half-lives using the pn-QRPA theory pi. The 
Q-values used in the calculations of rates in columns 3 and 4 are from [^9| while the mass formula of [[30| 
was used to calculate the rates in column 5. 55 Zn is an unstable nucleus according to the calculations of 
. All half-lives are given in seconds. Exponents are given in parenthesis. 



^4 


Element 


T c f l (terrestrial) 

( i Trainee tttiyti POI 


T°f l (stellar) 

I ) iralnpc TT"r~\m IvQI 
VdlLltJo 11 UI11 1— -'I 


T°f l f stellar 1 

( | "\7"CllnAtt TT"f~\m 1^(1 
V.^ VcxlLlco 11 Ulll 10L|I 


A O 




1.10 (-L ) 


i.iy (-L ) 


l.oo (-L) 






L.Li [-LJ 


L.QL y-L J 


L.yi \-L) 


45 


Fe 


7.03 (-3) 


9.37 (-3) 


1.01 (-2) 


4G 


Fe 


1.53 (-2) 


4.53 (-2) 


5.79 (-2) 


47 


Fe 


2.23 (-2) 


2.12 (-2) 


2.28 (-2) 


50 


Ni 


1.45 (-2) 


7.44 (-3) 


1.59 (-2) 


51 


Ni 


1.82 (-2) 


1.47 (-2) 


2.45 (-2) 


54 


Zn 


4.39 (-3) 


2.58 (-3) 


4.98 (-3) 


55 


Zn 


8.16 (-3) 


6.70 (-3) 




58 


Ge 


1.12 (-2) 


1.18 (-2) 


1.36 (-2) 


09 


Kr 


1.65 (-2) 


6.74 (-2) 


9.37 (-2) 


73 


Sr 


2.26 (-2) 


5.32 (-2) 


7.34 (-2) 


7G 


Sr 


1.22 (+1) 


4.91 (+1) 


4.68 (+1) 


80 


Zr 


4.72 (+0) 


7.24 (+0) 


4.64 (+0) 


84 


Mo 


3.07 (+0) 


3.40 (+0) 


2.85 (+0) 


87 


Ru 


3.50 (-2) 


6.25 (-2) 


8.14 (-2) 


91 


Pd 


3.30 (-2) 


4.51 (-2) 


8.69 (-2) 


92 


Pd 


1.24 (+0) 


4.32 (+0) 


9.79 (+0) 


93 


Pd 


2.51 (-1) 


3.15 (-1) 


2.53 (-1) 


95 


Ag 


9.62 (-1) 


8.31 (-1) 


9.86 (-1) 


9G 


Cd 


2.53 (-1) 


6.19 (-1) 


1.00 (+0) 


99 


In 


3.25 (+0) 


7.67 (+0) 


3.92 (+0) 
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Comparison of the terrestrial and stellar half-lives at lowest temperature and density 

Table (5): Comparison of calculated j3~ decay half-lives at temperature (T9 = 0.01) and density (pY e 
= 10 gem -3 ), with the calculation of terrestrial half-lives using the pn-QRPA theory [||. The Q- values 
used in the calculations of rates in columns 3 and 4 are from |29j while the mass formula of |5(| was 
used to calculate the rates in column 5. All half-lives are given in seconds. Exponents arc given in 
parenthesis. 



A Element T{f 2 (terrestrial) T{f % (stellar) T{f 2 (stellar) 

Q-values from [^9| Q- values from j29| Q-values from |3(J 



60 


Ti 


2.93 


(-3) 


3.43 


(-3) 


7.26 


(-4) 


Gl 


Ti 


1.43 


(-3) 


2.35 


(-3) 


2.64 


(-3) 


62 


V 


1.89 


(-2) 


1.94 


(-2) 


1.03 


(-1) 


63 


V 


4.99 


(-2) 


2.16 


(-2) 


3.47 


(-2) 


64 


Cr 


1.28 


(-2) 


5.02 


(-2) 


5.09 


(-3) 


69 


Pe 


4.97 


(-2) 


9.48 


(-2) 


3.51 ( 


+0) 


70 


Co 


4.05 


(-2) 


1.78 


(-2) 


1.45 


(-2) 


71 


Co 


4.26 


(-2) 


2.23 


(-2) 


2.02 


(-2) 


72 


Co 


2.23 


(-2) 


1.77 


(-2) 


1.70 


(-1) 


73 


Ni 


6.27 


(-1) 


2.39 


(-1) 


3.95 


(-1) 


74 


Ni 


5.63 


(-1) 


4.31 


(-1) 


1.21 ( 


+0) 


75 


Ni 


2.11 


(-1) 


1.81 


(-1) 


2.92 


(-1) 


76 


Ni 


1.17 


(-1) 


1.24 


(-1) 


2.12 


(-1) 


77 


Ni 


8.27 


(-2) 


9.77 


(-2) 


1.41 


(-1) 


78 


Ni 


3.96 


(-2) 


5.93 


(-2) 


7.60 


(-2) 


82 


Zn 


4.45 


(-2) 


8.97 


(-2) 


8.09 


(-2) 


85 


Ge 


1.31 


(-1) 


1.67 


(-1) 


1.30 


(-1) 


86 


Ge 


5.70 


(-2) 


7.04 


(-2) 


5.31 


(-2) 


88 


As 


5.70 


(-2) 


8.14 


(-2) 


6.53 


(-2) 


89 


As 


4.78 


(-2) 


5.08 


(-2) 


4.84 


(-2) 


96 


Kr 


2.67 


(-1) 


4.08 


(-1) 


3.34 


(-1) 


97 


Kr 


1.70 


(-1) 


3.29 


(-1) 


2.27 


(-1) 
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Accuracy of the pn-QRPA theory 



Table (6): The accuracy of the pn-QRPA model compared to experimental data (/3 + /EC decay). N 
denotes the number of experimentally known half-lives shorter than the limit in the second column, n is 
the number (and percentage) of isotopes reproduced under the condition given in the first column, and 
x is the average deviation defined in the text (see B for details). 



Conditions 


tt;» < 


N 


n 


n(%) 


X 


\/xi < 10 


10 e 


894 


706 


79.0 


2.057 




GO 


327 


304 


93.0 


1.718 




1 


81 


78 


96.3 


1.848 


Vie, < 2 


10 b 
GO 
1 


894 
327 
81 


489 
245 
59 


54.7 
74.9 
72.8 


1.363 
1.308 
1.230 



Table (7): The accuracy of the pn-QRPA model compared to experimental data (j3~ decay). N 
denotes the number of experimentally known half-lives shorter than the limit in the second column, n is 
the number (and percentage) of isotopes reproduced under the condition given in the first column, and 
x is the average deviation defined in the text (see for details) . 



Conditions 


T$(s) < 


N 


n 


n(%) 


X 


Vz; < 10 


10 B 
60 
1 


654 
325 
106 


472 
313 
105 


72.2 
96.3 
99.1 


1.85 ± 1.21 
1.67 ± 1.02 
1.44 ± 0.40 


Vxi < 5 


10 e 
60 
1 


654 
325 
106 


456 
307 
105 


69.7 
94.5 
99.1 


1.68 ± 0.76 
1.56 ± 0.66 
1.44 ± 0.40 


Vxi < 3 


10 e 
GO 
1 


654 
325 
106 


420 
295 
105 


64.2 
90.8 
99.1 


1.50 ± 0.46 
1.46 ± 0.43 
1.44 ± 0.40 


Vzi < 2 


10 e 
60 
1 


654 
325 
106 


369 
267 
96 


56.4 
82.2 
90.6 


1.37 ± 0.29 
1.36 ± 0.29 
1.35 ± 0.27 
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Predictive power of the pn-QRPA theory 



Table (8): Comparison of calculated (3 + / {EC) decay half-lives using the pn-QRPA theory and Q- values 
from H] , with the measured half- lives reported after the calculations done in [|| . All half- lives are given 
in seconds. Experimental half-lives are taken from [43-75] 



A 


Element 


rpexp 
1/2 


rpcal 
± l/2 


rpexp Irpcal 
1/2/1/2 


Comments 


70. 

t y 


v 


i a son -i- f{ nnn 

14.oUU zc U.UUU 


1 n 9fi/L 
lU.ZU^l 


u.oy 




OU 


Mr. 
IVIO 


i q finn + 1 inn 
ly.uuu zn ±.iuu 


7fi A^ 


o.yu 




88 
oo 


Tr- 
1 C 


a Ann -\- n son 


1 . i OO 


n 97 
u.z / 




yu 


nu 


nnn -i- 1 nnn 
y.uuu m i.uuu 


90 1 99 
zy . ±zz 


3 9zt 
o.z^ 




■:) 1 


rvu 


7 ^nn -i- n snn 

( .OUU ZC U.ouu 




n 79 
u. / z 




1 on 


Tn 
111 


a i nn -i- n onn 

u.iuu in u.yuu 


Z. lOO 


u.oo 




1 OA 
1U4 


DD 


n nod -u n i sn 

U.OZU ZIZ U.loU 




9 8fi 
Z.oU 




1UO 


o □ 


i ^nn -i- n 9nn 

l.OUU ZIZ U.ZUU 




1 9Q 
i.zy 




iuy 


±C 


/i Rnn -i- n ^nn 

^.UUU ZIZ U.OUU 


9 Q^7 
z.yo * 


U.O^l 




1 1 n 

11U 


Tn 


i & finn -i- n Rnn 

lO.UUU ZIZ u.ouu 


98 fino 
zo.uuy 


1 ^zt 




1 08 

1UO 




n n^fi -i- n nnfi 

U.UOO ZlZ u.uuu 


n zin9 


111 
ii.i 


partly alpna emission 






n a^Ci -u n n9n 

U.UOU ZlZ u.uzu 


n soft 
u.oyu 


1 "?& 
l.oo 




111 
in 




o c:nn -i- n 9nn 

Z.OUU ZlZ u.zuu 


9 Q77 
z.y ( * 


1 1 Q 
i . iy 




119 
1 Iz 




o /ion -i- n 1 1 n 

0.4ZU zt U.11U 


k n^^ 
o.uuo 


1 /IS 




1 1 ^ 
1 10 




a ar\n _i_ n 9nn 
U.UUU ZIZ u.zuu 


infill 


1 fil 

l.Ul 




111 
111 


-A.e 


n 7/in -i- n 9nn 

U. * ^±U ZlZ U.ZUU 


U.Ulo 




partly alpna emission 


119 
1 IZ 


Y 

^ve 


o 7nn -i- n snn 

Z. ( UU ZlZ U.oUU 


^ n^9 
o.uoz 


1 15 

1 . lo 




1 1 1 
1 lo 


Yo 


9 7/in -i- n nsn 

Z. ( ^±U ZlZ U.UoU 


9 sno 
z.ouy 


i n 1 ? 

l.UO 




118 
1 lo 


Dd 


c; 9nn -i- n 9nn 

O.ZUU ZlZ u.zuu 


Q AA7 


1 89 
l.oZ 






Ft 
£1/1 


i 7nn -i- n finn 

1. / UU ZIZ u.uuu 


^ n^9 


1 78 
1. ( o 




1 




n n^^ -i- n m n 

U.UOU ZIZ U.U1U 


n 

U. ±Oo 


o.yu 


pai Liy pioLoii emission 


1 £1 
101 


W 


U.41U zt U.U4U 




o.ou 


partly alpha emission 


162 


W 


1.390 ± 0.400 


4.325 


3.11 


partly alpha emission 


163 


w 


2.750 ± 0.250 


3.897 


1.42 


partly alpha emission 


164 


w 


6.400 ± 0.800 


10.273 


1.61 




165 


w 


5.100 ± 0.500 


8.748 


1.72 




166 


w 


18.800 ± 0.400 


21.117 


1.12 




167 


w 


19.900 ± 0.500 


24.663 


1.24 


partly alpha emission 


162 


Re 


0.100 ± 0.300 


0.514 


5.14 


partly alpha emission 


163 


Re 


0.260 ± 0.040 


1.132 


4.35 


partly alpha emission 


164 


Re 


0.880 ± 0.240 


0.867 


0.99 


partly alpha emission 


165 


Re 


2.400 ± 0.600 


2.018 


0.84 




166 


Re 


2.800 ± 0.300 


1.998 


0.71 


partly alpha emission 


167 


Re 


3.400 ± 0.400 


3.693 


1.09 


partly alpha emission 


168 


Re 


4.400 ± 0.100 


3.922 


0.89 




169 


Re 


8.100 ± 0.500 


10.369 


1.28 




169 


Ir 


0.400 ± 0.100 


1.523 


3.81 


partly alpha emission 


170 


Ir 


1.050 ± 0.150 


0.417 


0.40 


partly alpha emission 


171 


Ir 


1.500 ± 0.100 


2.436 


1.62 


partly alpha emission 


172 


Ir 


4.400 ± 0.300 


2.903 


0.66 
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Predictive power of the pn-QRPA theory 

Table (8) (contd.): 



A 


Element 


rpexp 


rpcal 


rpexp Irpcal 
1/2 / 1/2 


Comments 


174 


Ir 


9.000 ± 2.000 


5.546 


0.62 






175 


Ir 


9.000 ± 2.000 


10.652 


1.18 






176 


Ir 


8.000 ± 1.000 


16.729 


2.09 






177 


Ir 


30.000 ± 2.000 


42.039 


1.40 






174 


Au 


0.120 ± 0.020 


0.503 


4.19 


partly alpha 


emission 


175 


Au 


0.200 ± 0.022 


1.079 


5.39 


partly alpha 


emission 


176 


Au 


1.250 ± 0.300 


0.881 


0.70 


partly alpha 


emission 


177 


Au 


1.180 ± 0.070 


2.302 


1.95 


partly alpha 


emission 


178 


Au 


2.600 ± 0.500 


2.621 


1.01 


partly alpha 


emission 


179 


Au 


7.100 ± 0.300 


5.732 


0.81 


TtilT"! \\T C A XWn'A 




180 


Tl 


1.900 ± 0.900 


0.742 


0.39 


partly alpha 


emission 


181 


Tl 


3.400 ± 0.600 


2.331 


0.69 


partly alpha 


emission 


182 


Tl 


3.100 ± 1.000 


1.810 


0.58 






198 


At 


4.200 ± 0.300 


5.172 


1.23 


partly alpha 


emission 


234 


Am 


139.200 ± 4.800 


163.325 


1.17 






238 


Bk 


144.000 ± 5.000 


85.205 


0.59 






240 


Bk 


288.000 ±48.000 


357.840 


1.24 






239 


Cf 


39.000 ±37.000 


45.802 


1.17 


partly alpha 


emission 


240 


Cf 


63.600 ± 9.000 


56.025 


0.88 


partly alpha 


emission 


241 


Cf 


226.800 ±42.000 


184.563 


0.81 


partly alpha 


emission 


242 


Cf 


209.400 ± 7.200 


311.928 


1.49 






242 


Es 


7.000 ± 0.000 


15.271 


2.18 






249 


Fm 


156.000 ±42.000 


199.917 


1.28 






250 


Fm 


1800.00 ±180.00 


915.150 


0.51 


partly alpha 


emission 


253 


No 


102.000 ±18.000 


44.769 


0.44 


partly alpha 


emission 


254 


Lr 


13.000 ± 2.000 


24.733 


1.90 


partly alpha 


emission 


258 


Ha 


20.000 ±10.000 


13.857 


0.69 
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Predictive power of the pn-QRPA for proton-rich nuclei using different mass formulae 

Table (9): Comparison of calculated f3 + /(EC) decay half-lives using the pn-QRPA theory, with the 
measured half-lives reported after the calculations done in M. Columns 4-6 represent the calculated 
half-lives using the mass formulae of J7(|, and p9| ], respectively. The last three columns show the 
ratio of calculated half-lives (using mass formulae from |76]], |77| and ^9|, respectively) to measured 
ones. All half-lives are given in seconds. Experimental half-lives are taken from [43-75] 



A 


Element 


rpCXp 
1/2 


rpcal 
1/2 

1^X1111 ) 


rpcal 
1/2 

I'firnnto'l 


rpcal 
1/2 

I IVTUllCl i 


rpcal IrpCXp 

1/2/ 1/2 

fFTilFt 
^niii j 


rpcal irpexp 
1/2/ 1/2 

( rirnnfpi 


rpcal irpexp 
1/2/ 1/2 


99 




u.uuu in u.uuo 


n nee. 


U.UOl 


u.uuo 


IV. o 




1 1 9 




rv 

L 1 


n noi -U n nn/i 
U.UZ1 zt U.UU4 


U.UO / 


U.U4 ( 


U.UZo 


9 7Q 


9 9A 


1 HQ 


44 


Cr 


0.053 ± 0.004 


0.083 


0.068 


0.041 


1.57 


1.28 


0.77 


46 


Mn 


0.041 ± 0.007 


0.021 


0.114 


0.000 


0.50 


2.79 


0.00 


46 


Fe 


0.020 ± 0.020 


0.031 


0.024 


0.015 


1.53 


1.21 


0.77 


47 


Fe 


0.027 ± 0.032 


0.057 


0.050 


0.022 


2.12 


1.87 


0.83 


50 


Fe 


0.150 ± 0.030 


0.125 


0.214 


0.000 


0.83 


1.43 


0.00 


52 


Ni 


0.038 ± 0.005 


0.026 


0.073 


0.000 


0.67 


1.93 


0.00 


61 


Ga 


0.150 ± 0.030 


0.227 


0.298 


0.000 


1.51 


1.99 


0.00 


61 


Ge 


0.040 ± 0.015 


0.056 


0.054 


0.044 


1.39 


1.36 


1.11 


62 


Ge 


0.110 ± 0.060 


0.086 


0.082 


0.083 


0.78 


0.75 


0.76 


94 


Ag 


0.010 ± 0.000 


0.023 


0.023 


0.020 


2.33 


2.30 


2.01 


95 


Ag 


2.000 ± 0.100 


0.842 


0.654 


0.962 


0.42 


0.33 


0.48 


100 


Sn 


1.000 ± 0.800 


2.018 


1.151 


0.882 


2.02 


1.15 


0.88 


101 


Sn 


3.000 ± 1.000 


9.765 


6.049 


4.417 


3.26 


2.02 


1.47 


114 


Ba 


0.400 ± 0.300 


0.457 


0.510 


0.292 


1.14 


1.27 


0.73 


115 


Ba 


0.400 ± 0.200 


0.835 


0.675 


0.424 


2.09 


1.69 


1.06 


116 


Ba 


1.350 ± 0.150 


1.613 


1.649 


1.033 


1.19 


1.22 


0.76 


122 


Ce 


8.700 ± 0.700 


3.660 


2.934 


2.101 


0.42 


0.34 


0.24 


121 


Pr 


1.400 ± 0.800 


0.784 


0.813 


0.560 


0.56 


0.58 


0.40 


127 


Pr 


15.100 ± 0.000 


6.291 


15.514 


9.269 


0.42 


1.03 


0.61 


132 


Sm 


4.000 ± 0.300 


2.736 


2.533 


1.477 


0.68 


0.63 


0.37 
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Predictive power of the pn-QRPA for neutron-rich nuclei using different mass formulae 

Table (10): Comparison of calculated f3~ decay half-lives using the pn-QRPA theory, with the measured 
half-lives reported after the calculations done in S . Columns 4-6 represent the calculated half-lives using 
the mass formulae of J7(| , Q and [ p9|, r espectively. The last three columns show the ratio of calculated 
half-lives (using mass formulae from |76|, ]77| ] and ]2£|, respectively) to measured ones. All half- lives are 
given in seconds. Experimental half-lives are taken from [43-75] 



A 


Element 


rpCXp 
1/2 


rpcal 
1/2 


rpcal 
1/2 

I VJl UUtt I 


rpcal 
1/2 

I 1V_LU111_.I 1 


rpcal irpexp 
1/2/ 1/2 


rpcal IrpCXp 

1/2/ 1/2 

1 VJlUUlL 1 


rpcal irpexp 
1/2/ 1/2 

("Mollerl 


91 


NT 


n nor: _i_ pi m o 
u.uyo zz u.uio 


n 901 


901 


000 

U.UUU 


9 1 9 


9 1 9 


00 
u.uu 


99 


>J 


\J.\J£ L ± HZ U.UUI 




0^0 


000 
u.uuu 


9 D7 


9 OS 
. uo 


00 
u.uu 






U.UOZi ZZ U.UOI 


D 09^ 


09^ 


01 ^ 


n 98 


98 

U. ZiO 


Ifi 

U . J-U 


OA 




n Ofil 4- 09fi 


n m 9 


01 1 


007 
u.uu / 


n 9n 


1 7 
U. 1 / 


19 


97 


Np 
in c 


U.UJZi ZZ U.UUZi 




04fi 


09^ 

U. UZnJ 


1 


1 4^ 

_1_ .4:0 


7Q 

U. ( V 


98 
— o 


Np 

-N (. 


014 4- 01 

U.UJ-4: ZZ U.U1U 




09fi 


0^9 

U. UOZi 


9 A"\ 


1 8^ 
_L . OO 


9 9fi 

_ . _ U 


4'-! 
-to 


p 


U.Udd ZZ U.UUu 


\J.\JOO 


U . ± ± U 


034 


9 


3 48 

O .4:0 


1 09 


49 


o 


U . UUU ZZ u.uuu 




4 ^30 

4: . UOU 


1 ^0 


7 Qfi 


& OQ 


9 41 

Zi .4:± 




o 


n n&9 4- o 01 ^ 

u.uoz ZZ u.uio 


9QQ 


^Q7 
u.oy * 


1 0^ 

U. 1UO 


o.uo 


4 84 


1 98 


••1 A 
44 


Kyi 


U.4oU zt u.uou 


q i on 


o.UZU 


U.DDU 


7/19 
( .4/ 


11 .0 f 


1 .00 


45 


CI 


0.400 ± 0.043 


1.670 


1.670 


0.391 


4.17 


4.17 


0.98 


4G 


CI 


0.220 ± 0.040 


0.818 


0.848 


0.174 


3.72 


3.85 


0.79 


G5 


Fe 


0.450 ± 0.150 


0.335 


0.335 


0.210 


0.74 


0.74 


0.47 


68 


Fe 


0.100 ± 0.060 


0.160 


0.152 


0.209 


1.60 


1.52 


2.09 


68 


Co 


0.180 ± 0.100 


0.289 


0.289 


0.116 


1.61 


1.61 


0.64 


69 


Co 


0.270 ± 0.050 


0.071 


0.074 


0.082 


0.26 


0.27 


0.30 


71 


Co 


0.200 ± 0.050 


0.051 


0.059 


0.043 


0.25 


0.29 


0.21 


71 


Ni 


1.860 ± 0.350 


5.440 


5.860 


2.850 


2.92 


3.15 


1.53 


72 


Ni 


2.060 ± 0.300 


1.560 


1.680 


3.450 


0.76 


0.82 


1.67 


73 


Ni 


0.900 ± 0.150 


0.806 


1.040 


0.627 


0.90 


1.16 


0.70 


74 


Ni 


1.100 ± 0.500 


0.462 


0.543 


0.563 


0.42 


0.49 


0.51 


77 


Cu 


0.469 ± 0.008 


0.730 


1.080 


0.448 


1.56 


2.30 


0.96 


78 


Cu 


0.342 ± 0.008 


0.361 


0.653 


0.261 


1.06 


1.91 


0.76 


79 


Cu 


0.188 ± 0.011 


0.390 


0.608 


0.284 


2.07 


3.23 


1.51 


81 


Zn 


0.290 ± 0.050 


0.336 


0.665 


0.055 


1.16 


2.29 


0.19 


84 


Ga 


0.085 ± 0.010 


0.194 


0.406 


0.143 


2.28 


4.78 


1.68 


85 


Ge 


0.540 ± 0.050 


0.355 


0.761 


0.131 


0.66 


1.41 


0.24 


103 


Y 


0.230 ± 0.030 


0.162 


0.149 


0.362 


0.70 


0.65 


1.57 


105 


Zr 


1.000 ± 0.000 


0.820 


0.654 


0.241 


0.82 


0.65 


0.24 


107 


Nb 


0.330 ± 0.050 


0.602 


0.971 


0.271 


1.82 


2.94 


0.82 
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Predictive power of the pn-QRPA for neutron-rich nuclei using different mass formulae 
Table (10) (contd.): 



A 


Element 


rpexp 
1/2 


rpcal 
1/2 

("Hilf") 


rpcal 
1/2 

(Grootel 


rpcal 
1/2 

fMollorl 


rpcal JrpCXp 

1/2/ 1/2 

(HUH 


rpcal irpexp 
1/2/ 1/2 


rpcal irpexp 
1/2/ 1/2 

(Moller^ 


1 08 


Nh 


1 70 + 090 


389 


688 


1 87 


2.25 


4 05 


1 10 


109 


Nb 


n i on -)- n n^O 


0.542 


0.520 


0.129 


2.85 


2.74 


0.68 


1 1 


Nh 


i 7fi _|_ n non 

U.1IU —1- U.U^U 


984 

\J . .£04: 


116 


097 


1 67 


68 


57 

yj . O i 


1 0Q 


Mo 


n wn + n OfiO 


407 


809 


379 


77 


1 51 


79 


110 


Mo 




0.675 


0.593 


0.276 


2.25 


1.98 


0.92 


113 


Tc 


1 30 + 050 


0.244 


0.281 


0.095 


1.88 


2.16 


0.73 


115 


Ru 




0.434 


0.443 


0.240 


1.08 


1.11 


0.60 


1 1 7 


Rh 


n 44f) + 040 


408 


910 


0.226 


93 


48 


51 


1 1 Q 


Pd 


n qon 4 1 30 


1 900 


3 160 


577 


1 30 


3 43 


63 


120 


Pd 


0.500 ± 0.100 


0.695 


1.660 


0.232 


1.39 


3.32 


0.46 


125 


Ag 


0.156 ± 0.007 


0.164 


0.249 


0.117 


1.05 


1.60 


0.75 


126 


Ag 


0.097 ± 0.008 


0.140 


0.194 


0.090 


1.44 


2.00 


0.93 


127 


Ag 


0.109 ± 0.015 


0.092 


0.151 


0.089 


0.84 


1.39 


0.81 


149 


Ba 


0.344 ± 0.007 


0.287 


0.345 


0.365 


0.83 


1.00 


1.06 


150 


La 


0.860 ± 0.050 


0.771 


1.460 


0.405 


0.90 


1.70 


0.47 


185 


Hf 


210.000 ±36.000 


584.000 


766.000 


204.000 


2.78 


3.65 


0.97 


205 


Au 


31.000 ± 2.000 


35.700 


50.000 


25.100 


1.15 


1.61 


0.81 


216 


Bi 


216.000 ±24.000 


189.000 


191.000 


223.000 


0.88 


0.88 


1.03 


221 


At 


138.000 ± 12.000 


256.000 


462.000 


564.000 


1.86 


3.35 


4.09 


222 


At 


54.000 ± 10.000 


63.300 


155.000 


10.800 


1.17 


2.87 


0.20 


223 


At 


50.000 ± 7.000 


41.600 


84.100 


70.700 


0.83 


1.68 


1.41 


232 


Fr 


5.000 ± 1.000 


1.690 


15.900 


0.903 


0.34 


3.18 


0.18 


233 


Ra 


30.000 ± 5.000 


86.000 


622.000 


52.000 


2.87 


20.73 


1.73 


234 


Ra 


30.000 ±10.000 


34.300 


127.000 


42.000 


1.14 


4.23 


1.40 


237 


Th 


300.000 ± 54.000 


221.000 


389.000 


311.000 


0.74 


1.30 


1.04 
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Table (11): Comparison of calculated B(GT) + strength function using the pn-QRPA theory (this 
work), with the measured ones from (n,p) reactions. QRPA and SM represent calculations using the 
pn-QRPA theory and those from the recent large-scale shell model calculations |7i| , respectively. The 
experimental data are taken from [79-83] . 



Nucleus 


Measured Strength 


QRPA 


SM 


Sly 


1.2 ± 0.1 


1.4 


1.4 


54 Fe 


3.5 ± 0.3 


3.8 


3.6 


55 Mn 


1.7 ± 0.2 


2.4 


2.2 


56 Fe 


2.9 ± 0.3 


2.7 


2.7 


58 Ni 


3.8 ± 0.4 


3.3 


4.4 


59 Co 


1.9 ± 0.1 


0.8 


2.5 


60 Ni 


3.1 ± 0.1 


7.0 


3.4 


62 Ni 


2.5 ± 0.1 


4.4 


2.1 




1.7 ± 0.1 


3.G 


1.3 


70 Ge 


0.7 ± 0.1 


0.3 




72 Ge 


0.2 ± 0.1 


0.1 





Table (12): Comparison of calculated B(GT)~ strength function using the pn-QRPA theory (this 
work), with the measured ones from (p,n) reactions. QRPA and SM represent calculations using the 
pn-QRPA theory and those from the recent large-scale shell model calculations |7§], respectively. The 
experimental data are taken from [84-87]. See |84j] regarding the uncertainty in measured Gamow- Teller 
strength value of 32 S. 



Nucleus 


Measured Strength 


QRPA 


SM 




2.21 


1.1 




Sly 


12.6 ± 2.5 


16.6 




54 Fe 


7.8 ± 1.9 


5.9 


6.9 


58 Ni 


7.4 ± 1.8 


7.5 


7.7 


eo Ni 


7.2 ± 1.8 


7.7 


10.0 


71 Ga 


4.3 ± 0.7 


5.6 
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EXPLANATION OF TABLE 



TABLE A. fp-Shell Nuclei Weak Rates in Stellar Matter 



The calculated weak interaction rates [Eqs. (16) and (19)-(21)] are all tabulated in logioA. The 
probabilities of /3-delayed proton (neutron) emission [Eq. (22)] are also tabulated in logarithmic scale. 
These probabilities are calculated only to one significant figure and are given up to three places of 
decimal only for designing purposes. All rates listed for a particular direction concern the parent 
nucleus except for the last two columns which concern the daughter nucleus. For each daughter nucleus, 
cither the proton energy rate and probability of /3-delayed proton emission is stated (if S p < S n ) or the 
neutron energy rate and probability of /3-delayed neutron emission is stated (if S n < S p ). In the table 
-100 means that the rate (or the probability) is smaller than 10~ 100 . 



Q 


Mass of parent minus mass of daughter nucleus 




SnP 


Separation energy of neutron for parent nucleus 




SpP 


Separation energy of proton for parent nucleus 




SnD 


Separation energy of neutron for daughter nucleus 




SpD 


Separation energy of neutron for daughter nucleus 




ADen 


log (pY e ) (g.cm~ 3 ), where p is the density of the baryon, and Y e 
electron number to the baryon number 


is the ratio of the 


T9 


Temperature in units of 10 9 K 




EFermi 


Total Fermi energy of electron and positron, including the rest mass (McV) 


E+Cap 


Positron capture rate (s _1 ) 
Electron decay rate (s _1 ) 




E-Dec 




ANuEn 


Anti-Neutrino energy loss rate (MeV.s -1 ) 

Gamma ray heating rate (MeV.s -1 ), tabulated separately for /3 + 
and /3~ direction 




GamEn 


direction 


E-Cap 


Electron capture rate (s _1 ) 
Positron decay rate (s _1 ) 




E+Dec 




NuEn 


Neutrino energy loss rate (MeV.s -1 ) 




ProEn 


Energy rate of /3-delayed proton (MeV.s -1 ) 




NeuEn 


Energy rate of /3-delayed neutron (MeV.s -1 ) 




PPEm 


Probability of /3-delayed proton emission 




PNEm 


Probability of /3-delayed neutron emission 
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40 
20 



CA 



40 cr 

21 



Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-76.256 


-100 


-76.205 


-80.721 


-85.961 


-9.000 


3.0 


3.0 


-25.504 


-100 


-25.348 


-26.538 


-28.464 


-3.000 


3.0 


30.0 


0.278 


-100 


1.277 


-0.468 


0.766 


-0.100 


7.0 


1.0 


-82.076 


-100 


-82.025 


-86.541 


-91.781 


-9.000 


7.0 


3.0 


-27.220 


-100 


-27.064 


-28.254 


-30.180 


-3.000 


7.0 


30.0 


0.275 


-100 


1.274 


-0.471 


0.763 


-0.100 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-65.697 


-100 


-65.541 


-66.731 


-68.657 


-3.000 


11.0 


30.0 


-3.588 


-100 


-2.589 


-4.334 


-3.100 


-0.100 



40 
20 



CA 



20 



40 
19 



K 



21 



Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-25.483 


-25.143 


-25.002 


-25.103 


-56.523 


-30.000 


3.0 


3.0 


-9.824 


-9.645 


-9.280 


-9.484 


-19.330 


-9.000 


3.0 


30.0 


1.905 


-2.098 


3.004 


2.485 


1.260 


-0.800 


7.0 


1.0 


-19.663 


-25.143 


-20.234 


-19.283 


-50.703 


-30.000 


7.0 


3.0 


-8.162 


-9.640 


-8.021 


-7.840 


-17.615 


-9.000 


7.0 


30.0 


1.909 


-2.097 


3.007 


2.489 


1.263 


-0.800 


11.0 


1.0 


3.650 


-25.143 


4.781 


4.182 


2.745 


-0.700 


11.0 


3.0 


3.651 


-9.640 


4.783 


4.183 


2.747 


-0.700 


11.0 


30.0 


4.793 


-1.905 


6.112 


5.373 


4.458 


-0.600 








40 
21 


SC 19 — ► 


40 Tl 
22 11 18 







Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-81.245 


-100 


-81.194 


-100 


-81.298 


0.000 


3.0 


3.0 


-26.642 


-100 


-26.490 


-100 


-26.661 


0.000 


3.0 


30.0 


0.015 


-100 


0.994 


-100 


0.390 


0.000 


7.0 


1.0 


-87.065 


-100 


-87.014 


-100 


-87.118 


0.000 


7.0 


3.0 


-28.357 


-100 


-28.206 


-100 


-28.377 


0.000 


7.0 


30.0 


0.011 


-100 


0.991 


-100 


0.387 


0.000 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-66.834 


-100 


-66.683 


-100 


-66.854 


0.000 


11.0 


30.0 


-3.852 


-100 


-2.872 


-100 


-3.476 


0.000 








40 
21 


SC 19 — ► 


20 CA 20 







Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-3.236 


1.696 


2.374 


2.320 


-2.292 


-3.000 


3.0 


3.0 


-0.961 


1.664 


2.344 


2.307 


-1.816 


-3.000 


3.0 


30.0 


2.812 


1.673 


4.094 


3.516 


2.549 


-0.900 


7.0 


1.0 


0.474 


1.696 


2.427 


2.347 


-1.189 


-3.000 


7.0 


3.0 


0.455 


1.665 


2.397 


2.333 


-0.957 


-3.000 


7.0 


30.0 


2.816 


1.673 


4.097 


3.520 


2.552 


-0.900 


11.0 


1.0 


5.285 


1.696 


6.714 


5.895 


5.196 


-0.800 


11.0 


3.0 


5.262 


1.665 


6.690 


5.883 


5.221 


-0.800 


11.0 


30.0 


5.397 


1.746 


6.858 


6.030 


5.451 


-0.700 
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40 T | 
22 1 1 



Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


3.0 


3.0 


-36.419 


-100 


-36.269 


-100 


-36.314 


0.000 


3.0 


30.0 


-1.620 


-100 


-0.661 


-100 


-1.005 


0.000 


7.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


7.0 


3.0 


-38.135 


-100 


-37.985 


-100 


-38.029 


0.000 


7.0 


30.0 


-1.624 


-100 


-0.665 


-100 


-1.008 


0.000 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-76.612 


-100 


-76.462 


-100 


-76.507 


0.000 


11.0 


30.0 


-5.487 


-100 


-4.528 


-100 


-4.871 


0.000 



40 
22 



Tl 



18 



40 cr 
91 -J^- 



19 



Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-3.260 


1.452 


2.100 


1.257 


1.338 


-0.300 


3.0 


3.0 


-0.948 


1.452 


2.103 


1.257 


1.343 


-0.300 


3.0 


30.0 


2.851 


1.373 


4.087 


2.226 


3.420 


-0.090 


7.0 


1.0 


0.456 


1.452 


2.176 


1.279 


1.452 


-0.300 


7.0 


3.0 


0.473 


1.452 


2.181 


1.279 


1.459 


-0.300 


7.0 


30.0 


2.855 


1.373 


4.090 


2.229 


3.423 


-0.090 


11.0 


1.0 


5.435 


1.452 


6.827 


4.676 


6.087 


-0.060 


11.0 


3.0 


5.436 


1.452 


6.828 


4.676 


6.087 


-0.060 


11.0 


30.0 


5.484 


1.452 


6.915 


4.711 


6.149 


-0.060 










20 CA 21 


41 cr 
—> 21 ov - 


20 





Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-37.297 


-43.398 


-37.245 


-41.481 


-53.133 


-20.000 


3.0 


3.0 


-12.955 


-16.868 


-12.798 


-14.400 


-18.304 


-5.000 


3.0 


30.0 


0.894 


-5.663 


1.920 


-0.086 


1.667 


-0.080 


7.0 


1.0 


-43.116 


-45.778 


-43.065 


-45.784 


-58.953 


-20.000 


7.0 


3.0 


-14.671 


-17.505 


-14.513 


-16.100 


-20.019 


-5.000 


7.0 


30.0 


0.891 


-5.665 


1.916 


-0.090 


1.664 


-0.080 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-53.148 


-55.860 


-52.990 


-54.572 


-58.497 


-5.000 


11.0 


30.0 


-2.970 


-9.168 


-1.945 


-3.945 


-2.198 


-0.080 



41 
20 



CA 



21 



19 



K 



22 



Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-7.534 


-20.309 


-7.712 


-12.560 


-49.416 


-40.000 


3.0 


3.0 


-4.649 


-8.688 


-4.515 


-6.374 


-17.128 


-10.000 


3.0 


30.0 


1.867 


-1.706 


2.934 


2.325 


2.072 


-0.300 


7.0 


1.0 


-3.438 


-20.309 


-3.405 


-8.473 


-43.596 


-40.000 


7.0 


3.0 


-3.097 


-8.687 


-2.907 


-4.822 


-15.415 


-10.000 


7.0 


30.0 


1.870 


-1.705 


2.937 


2.328 


2.075 


-0.300 


11.0 


1.0 


4.374 


-20.309 


5.557 


5.034 


3.455 


-1.000 


11.0 


3.0 


4.378 


-8.687 


5.563 


5.033 


3.478 


-1.000 


11.0 


30.0 


4.873 


-1.569 


6.150 


5.221 


5.239 


-0.200 
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41 
21 



sc 



41 T i 

22 1 ' 



Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-69.473 


-100 


-69.423 


-75.817 


-93.550 


-20.000 


3.0 


3.0 


-23.545 


-100 


-23.393 


-26.513 


-30.873 


-8.000 


3.0 


30.0 


-0.755 


-100 


0.226 


-2.498 


-0.581 


-0.200 


7.0 


1.0 


-75.293 


-100 


-75.242 


-81.637 


-99.369 


-20.000 


7.0 


3.0 


-25.261 


-100 


-25.109 


-28.228 


-32.588 


-8.000 


7.0 


30.0 


-0.758 


-100 


0.223 


-2.502 


-0.584 


-0.200 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-63.738 


-100 


-63.586 


-66.706 


-71.066 


-8.000 


11.0 


30.0 


-4.621 


-100 


-3.640 


-6.364 


-4.447 


-0.200 



41 
21 



SC 



20 



41 
20 



CA 



21 



Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-4.416 


-0.112 


0.355 


-3.080 


-22.675 


-20.000 


3.0 


3.0 


-2.163 


-0.181 


0.296 


-1.415 


-9.162 


-8.000 


3.0 


30.0 


1.938 


-0.192 


3.116 


2.300 


1.416 


-0.900 


7.0 


1.0 


-0.691 


-0.112 


0.566 


-2.071 


-16.855 


-20.000 


7.0 


3.0 


-0.740 


-0.180 


0.512 


-1.241 


-7.447 


-7.000 


7.0 


30.0 


1.942 


-0.191 


3.120 


2.303 


1.419 


-0.900 


11.0 


1.0 


4.445 


-0.112 


5.786 


3.768 


4.316 


-0.500 


11.0 


3.0 


4.416 


-0.180 


5.748 


4.097 


4.329 


-0.500 


11.0 


30.0 


4.694 


-0.071 


6.065 


5.101 


4.578 


-0.600 










41 Tl 

22 11 19 


^ 23 v 


18 





Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-82.127 


-100 


-82.076 


-100 


-83.571 


0.000 


3.0 


3.0 


-27.435 


-100 


-27.283 


-100 


-28.464 


0.000 


3.0 


30.0 


-0.805 


-100 


0.179 


-100 


-0.401 


0.000 


7.0 


1.0 


-87.947 


-100 


-87.896 


-100 


-89.391 


0.000 


7.0 


3.0 


-29.150 


-100 


-28.998 


-100 


-30.180 


0.000 


7.0 


30.0 


-0.809 


-100 


0.176 


-100 


-0.404 


0.000 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-67.628 


-100 


-67.476 


-100 


-68.657 


0.000 


11.0 


30.0 


-4.672 


-100 


-3.687 


-100 


-4.267 


0.000 



41 
22 



Tl 



19 



41 Cf 
21 3 ^ 



20 



Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-3.591 


0.944 


1.581 


0.759 


1.241 


-0.300 


3.0 


3.0 


-1.275 


0.940 


1.581 


0.754 


1.244 


-0.300 


3.0 


30.0 


2.686 


0.972 


3.887 


1.321 


3.547 


-0.020 


7.0 


1.0 


0.132 


0.944 


1.681 


0.774 


1.362 


-0.300 


7.0 


3.0 


0.149 


0.941 


1.685 


0.769 


1.367 


-0.300 


7.0 


30.0 


2.689 


0.973 


3.891 


1.324 


3.551 


-0.020 


11.0 


1.0 


5.286 


0.944 


6.645 


3.959 


6.139 


-0.020 


11.0 


3.0 


5.282 


0.941 


6.642 


3.952 


6.137 


-0.020 


11.0 


30.0 


5.370 


1.054 


6.774 


3.756 


6.292 


-0.008 
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42 
20 



CA 



42 cr 

21 3 ^ 



Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-35.853 


-47.080 


-36.198 


-38.323 


-63.872 


-30.000 


3.0 


3.0 


-12.594 


-16.777 


-12.443 


-14.215 


-21.477 


-8.000 


3.0 


30.0 


1.266 


-1.940 


2.323 


1.319 


1.544 


-0.200 


7.0 


1.0 


-41.672 


-48.803 


-42.017 


-44.143 


-69.692 


-30.000 


7.0 


3.0 


-14.309 


-16.933 


-14.158 


-15.882 


-23.192 


-8.000 


7.0 


30.0 


1.263 


-1.941 


2.319 


1.315 


1.540 


-0.200 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-52.787 


-52.818 


-52.478 


-52.986 


-61.669 


-8.000 


11.0 


30.0 


-2.595 


-5.293 


-1.541 


-2.542 


-2.321 


-0.200 








42 
20 


CA 22 - 


42 w 
19 K 23 







Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-26.477 


-30.847 


-26.972 


-28.640 


-67.776 


-40.000 


3.0 


3.0 


-10.901 


-12.362 


-10.818 


-11.313 


-22.919 


-10.000 


3.0 


30.0 


1.482 


-1.062 


2.568 


2.032 


1.138 


-0.800 


7.0 


1.0 


-21.464 


-30.847 


-21.772 


-22.838 


-61.956 


-40.000 


7.0 


3.0 


-9.256 


-12.347 


-9.174 


-9.639 


-21.203 


-10.000 


7.0 


30.0 


1.485 


-1.061 


2.572 


2.035 


1.142 


-0.800 


11.0 


1.0 


3.832 


-30.847 


4.962 


4.320 


1.484 


-3.000 


11.0 


3.0 


3.831 


-12.346 


4.963 


4.319 


1.487 


-3.000 


11.0 


30.0 


4.481 


-0.955 


5.760 


5.082 


4.332 


-0.600 








42 
21 


SC 21 - 


42 Tl 
22 11 20 







Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-39.232 


-100 


-39.181 


-60.032 


-64.586 


-20.000 


3.0 


3.0 


-13.288 


-100 


-13.131 


-19.841 


-21.648 


-8.000 


3.0 


30.0 


0.307 


-100 


1.316 


-0.267 


0.664 


-0.200 


7.0 


1.0 


-45.052 


-100 


-45.001 


-65.852 


-70.406 


-20.000 


7.0 


3.0 


-15.003 


-100 


-14.846 


-21.556 


-23.364 


-8.000 


7.0 


30.0 


0.303 


-100 


1.313 


-0.270 


0.661 


-0.200 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-53.481 


-100 


-53.323 


-60.034 


-61.841 


-8.000 


11.0 


30.0 


-3.559 


-100 


-2.550 


-4.132 


-3.202 


-0.200 



42 
21 



SC 



21 



42 TA 

20 



22 



Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-4.288 


-0.008 


0.453 


-2.936 


-31.563 


-30.000 


3.0 


3.0 


-2.116 


-0.182 


0.290 


-2.155 


-11.505 


-10.000 


3.0 


30.0 


1.918 


-1.118 


3.024 


2.590 


1.845 


-0.500 


7.0 


1.0 


-0.561 


-0.008 


0.673 


-1.421 


-25.743 


-20.000 


7.0 


3.0 


-0.690 


-0.181 


0.519 


-1.078 


-9.790 


-9.000 


7.0 


30.0 


1.921 


-1.117 


3.027 


2.594 


1.848 


-0.500 


11.0 


1.0 


4.733 


-0.008 


6.045 


5.256 


4.270 


-0.800 


11.0 


3.0 


4.722 


-0.181 


6.009 


5.322 


4.413 


-0.700 


11.0 


30.0 


4.829 


-0.996 


6.136 


5.403 


5.027 


-0.300 



40 









42 Tl 

22 1 1 20 — * 


42 \/ 

oq V 1 9 






Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-92.376 


-100 


-92.325 


-100 


-92.590 


0.000 


3.0 


3.0 


-30.760 


-100 


-30.609 


-100 


-30.951 


0.000 


3.0 


30.0 


-0.720 


-100 


0.251 


-100 


-0.083 


0.000 


7.0 


1.0 


-98.196 


-100 


-98.145 


-100 


-98.410 


0.000 


7.0 


3.0 


-32.475 


-100 


-32.325 


-100 


-32.666 


0.000 


7.0 


30.0 


-0.723 


-100 


0.248 


-100 


-0.087 


0.000 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-70.953 


-100 


-70.802 


-100 


-71.144 


0.000 


11.0 


30.0 


-4.587 


-100 


-3.615 


-100 


-3.950 


0.000 








42 
22 


Tl 20 -> 


42 cr 

21 21 






Aden 


T9 


E-Cap 


E+Dec 


NuEn GamEn ProEn PPEm 


3.0 


1.0 


-3.720 


0.569 


1.043 


0.167 - 


5.834 -6.000 


3.0 


3.0 


-1.400 


0.568 


1.053 


0.172 - 


3.081 -4.000 


3.0 


30.0 


2.843 


1.217 


4.050 


2.780 


2.985 -0.400 


7.0 


1.0 


0.002 


0.569 


1.256 


0.290 - 


1.774 -2.000 


7.0 


3.0 


0.024 


0.569 


1.270 


0.295 - 


1.517 -2.000 


7.0 


30.0 


2.847 


1.218 


4.054 


2.783 


2.988 -0.400 


11.0 


1.0 


5.201 


0.569 


6.552 


4.670 


5.230 -0.400 


11.0 


3.0 


5.202 


0.569 


6.553 


4.671 


5.231 -0.400 


11.0 


30.0 


5.526 


1.298 


6.931 


5.362 


5.848 -0.300 








43 (~ A 

20 LA 23 - 


* 2I SC 22 




Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-17.683 


-25.442 


-17.639 


-22.136 


-47.013 


-30.000 


3.0 


3.0 


-7.685 


-11.028 


-7.510 


-8.827 


-16.248 


-8.000 


3.0 


30.0 


1.515 


-2.279 


2.559 


1.456 


2.085 


-0.100 


7.0 


1.0 


-23.502 


-27.651 


-23.459 


-27.870 


-52.831 


-30.000 


7.0 


3.0 


-9.400 


-11.091 


-9.215 


-10.444 


-17.963 


-8.000 


7.0 


30.0 


1.512 


-2.280 


2.556 


1.453 


2.082 


-0.100 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-47.877 


-46.559 


-46.689 


-46.796 


-56.441 


-10.000 


11.0 


30.0 


-2.348 


-5.632 


-1.305 


-2.405 


-1.779 


-0.100 


43 r A _^ 43 1/ 
20 LA 23 —* 19 ^ 24 



Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-17.961 


-26.918 


-18.499 


-20.815 


-67.252 


-50.000 


3.0 


3.0 


-8.399 


-11.348 


-8.431 


-8.644 


-22.313 


-10.000 


3.0 


30.0 


1.533 


-2.313 


2.569 


2.128 


1.279 


-0.400 


7.0 


1.0 


-12.181 


-26.918 


-12.704 


-14.995 


-61.432 


-50.000 


7.0 


3.0 


-6.689 


-11.345 


-6.720 


-6.928 


-20.597 


-10.000 


7.0 


30.0 


1.537 


-2.312 


2.572 


2.132 


1.283 


-0.400 


11.0 


1.0 


4.039 


-26.918 


5.181 


4.650 


2.678 


-2.000 


11.0 


3.0 


4.114 


-11.345 


5.260 


4.730 


2.769 


-2.000 


11.0 


30.0 


4.659 


-2.160 


5.899 


5.240 


4.523 


-0.300 



41 



43 
21 



SC 



43 
22 



Tl 



Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-39.083 


-100 


-39.031 


-42.382 


-67.354 


-30.000 


3.0 


3.0 


-13.552 


-100 


-13.395 


-15.379 


-22.324 


-8.000 


3.0 


30.0 


0.656 


-100 


1.676 


0.532 


0.968 


-0.200 


7.0 


1.0 


-44.903 


-100 


-44.851 


-48.202 


-73.174 


-30.000 


7.0 


3.0 


-15.268 


-100 


-15.111 


-17.094 


-24.040 


-8.000 


7.0 


30.0 


0.652 


-100 


1.672 


0.529 


0.964 


-0.200 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-53.745 


-100 


-53.588 


-55.572 


-62.517 


-8.000 


11.0 


30.0 


-3.209 


-100 


-2.190 


-3.332 


-2.898 


-0.200 








43 
21 


SC 22 - 


i CA 23 







Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-6.480 


-4.382 


-4.553 


-5.516 


-40.994 


-40.000 


3.0 


3.0 


-3.775 


-3.301 


-2.949 


-4.230 


-14.461 


-10.000 


3.0 


30.0 


1.680 


-2.050 


2.787 


2.042 


1.575 


-0.400 


7.0 


1.0 


-2.658 


-4.381 


-2.254 


-3.279 


-35.174 


-30.000 


7.0 


3.0 


-2.319 


-3.296 


-1.761 


-3.117 


-12.746 


-10.000 


7.0 


30.0 


1.683 


-2.049 


2.791 


2.046 


1.578 


-0.400 


11.0 


1.0 


3.822 


-4.381 


4.987 


4.489 


3.656 


-0.600 


11.0 


3.0 


3.886 


-3.296 


5.067 


4.458 


3.679 


-0.500 


11.0 


30.0 


4.573 


-1.877 


5.887 


4.873 


4.731 


-0.300 










43 Tl 
22 1 1 21 


43 \/ 
-» 23 V 20 







Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-60.958 


-100 


-60.907 


-62.694 


-70.918 


-9.000 


3.0 


3.0 


-20.423 


-100 


-20.270 


-21.521 


-25.563 


-5.000 


3.0 


30.0 


-0.072 


-100 


0.922 


-1.417 


0.590 


-0.100 


7.0 


1.0 


-66.778 


-100 


-66.727 


-68.514 


-76.738 


-9.000 


7.0 


3.0 


-22.139 


-100 


-21.985 


-23.236 


-27.278 


-5.000 


7.0 


30.0 


-0.075 


-100 


0.919 


-1.420 


0.587 


-0.100 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-60.616 


-100 


-60.463 


-61.713 


-65.756 


-5.000 


11.0 


30.0 


-3.938 


-100 


-2.944 


-5.283 


-3.276 


-0.100 



43 
22 



Tl 



21 



43 
21 



SC 



22 



Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-4.213 


0.076 


0.567 


-1.249 


-6.527 


-6.000 


3.0 


3.0 


-1.955 


-0.043 


0.454 


-0.970 


-3.640 


-3.000 


3.0 


30.0 


2.439 


0.110 


3.575 


2.397 


2.983 


-0.100 


7.0 


1.0 


-0.472 


0.076 


0.766 


-0.479 


-2.409 


-2.000 


7.0 


3.0 


-0.518 


-0.043 


0.669 


-0.358 


-2.089 


-2.000 


7.0 


30.0 


2.443 


0.111 


3.579 


2.400 


2.987 


-0.100 


11.0 


1.0 


5.046 


0.076 


6.345 


5.187 


5.315 


-0.300 


11.0 


3.0 


5.030 


-0.043 


6.322 


5.182 


5.335 


-0.200 


11.0 


30.0 


5.258 


0.203 


6.601 


5.059 


5.937 


-0.070 



42 









43 \/ 
23 V 20 


oil C R 1 Q 






Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-90.162 


-100 


-90.111 


-90.160 


-100 


-20.000 


3.0 


3.0 


on 0"70 

-30.872 


-100 


-30.721 


-30.870 


-38.039 


-8.000 


3.0 


30.0 


-1.803 


-100 


-0.838 


-2.240 


-1.210 


-0.200 


7.0 


1.0 


-95.982 


-100 


-95.931 


-95.980 


-100 


-20.000 


7.0 


3.0 


-32.587 


-100 


-32.436 


-32.586 


-39.754 


-8.000 


7.0 


30.0 


-1.806 


-100 


-0.841 


-2.243 


-1.213 


-0.200 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-71.065 


-100 


-70.914 


-71.063 


-78.231 


-8.000 


11.0 


30.0 


-5.669 


-100 


-4.704 


-6.106 


-5.076 


-0.200 








43 
23 


V 20 ^ 


H Tl 2i 






Aden 


T9 


E-Cap 


E+Dec 


NuEn GamEn ProEn PPEm 


3.0 


1.0 


-3.665 


1.186 


1.870 


1.183 - 


0.728 -2.000 


3.0 


3.0 


-1.334 


1.212 


1.899 


1.211 - 


0.663 -2.000 


3.0 


30.0 


2.530 


1.153 


3.751 


2.218 


2.907 -0.300 


7.0 


1.0 


0.043 


1.186 


1.927 


1.208 - 


0.195 -2.000 


7.0 


3.0 


0.084 


1.212 


1.957 


1.235 - 


0.127 -2.000 


7.0 


30.0 


2.533 


1.154 


3.755 


2.221 


2.910 -0.300 


11.0 


1.0 


5.134 


1.186 


6.487 


4.631 


5.715 -0.200 


11.0 


3.0 


5.161 


1.212 


6.515 


4.657 


5.742 -0.200 


11.0 


30.0 


5.230 


1.224 


6.626 


4.704 


5.839 -0.200 








2Q CA 24 — 


44 cr 
» 21 - 51 - 23 




Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-23.429 


-29.193 


-23.850 


-26.210 


-62.297 


-40.000 


3.0 


3.0 


-9.852 


-11.225 


-9.735 


-10.124 


-20.513 


-10.000 


3.0 


30.0 


1.603 


-0.854 


2.671 


2.047 


1.835 


-0.300 


7.0 


1.0 


-29.249 


-29.513 


-29.363 


-31.949 


-67.996 


-40.000 


7.0 


3.0 


-11.564 


-11.428 


-11.174 


-11.800 


-22.228 


-10.000 


7.0 


30.0 


1.600 


-0.855 


2.668 


2.043 


1.832 


-0.300 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-50.042 


-46.815 


-46.955 


-48.255 


-60.706 


-10.000 


11.0 


30.0 


-2.258 


-3.967 


-1.190 


-1.811 


-2.029 


-0.300 


44 CA ^ 44 \s 

20 LA 24 — > 19 l\ 25 


Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-36.662 


-48.869 


-37.075 


-39.520 


-76.691 


-40.000 


3.0 


3.0 


-13.810 


-17.353 


-13.758 


-14.377 


-25.676 


-10.000 


3.0 


30.0 


1.447 


-1.434 


2.522 


1.937 


1.135 


-0.700 


7.0 


1.0 


-31.221 


-48.869 


-31.703 


-33.776 


-70.871 


-40.000 


7.0 


3.0 


-12.179 


-17.347 


-12.102 


-12.666 


-23.961 


-10.000 


7.0 


30.0 


1.450 


-1.433 


2.525 


1.940 


1.138 


-0.700 


11.0 


1.0 


2.666 


-48.869 


3.674 


3.178 


1.969 


-0.800 


11.0 


3.0 


2.663 


-17.347 


3.673 


3.175 


1.968 


-0.800 


11.0 


30.0 


4.501 


-1.301 


5.761 


5.050 


4.332 


-0.600 



43 



44 
21 



SC 



44 jl 
22 1 1 



Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-11.109 


-31.751 


-11.059 


-25.595 


-54.963 


-40.000 


3.0 


3.0 


-7.599 


-12.841 


-7.358 


-9.348 


-18.270 


-10.000 


3.0 


30.0 


1.381 


-3.550 


2.420 


1.827 


1.722 


-0.200 


7.0 


1.0 


-16.928 


-32.988 


-16.879 


-31.371 


-60.783 


-40.000 


7.0 


3.0 


-9.308 


-13.151 


-9.068 


-11.048 


-19.985 


-10.000 


7.0 


30.0 


1.377 


-3.552 


2.417 


1.824 


1.718 


-0.200 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-47.785 


-49.846 


-47.544 


-49.070 


-58.463 


-10.000 


11.0 


30.0 


-2.483 


-6.970 


-1.444 


-2.035 


-2.144 


-0.200 








44 
21 


SC 23 - 


20 CA 24 







Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-7.015 


-4.558 


-4.629 


-4.491 


-49.294 


-40.000 


3.0 


3.0 


-4.740 


-4.845 


-4.253 


-4.247 


-17.073 


-10.000 


3.0 


30.0 


1.807 


-1.466 


2.874 


2.530 


1.880 


-0.400 


7.0 


1.0 


-3.190 


-4.557 


-2.765 


-2.980 


-43.474 


-40.000 


7.0 


3.0 


-3.243 


-4.836 


-2.838 


-2.866 


-15.358 


-10.000 


7.0 


30.0 


1.811 


-1.465 


2.878 


2.534 


1.883 


-0.400 


11.0 


1.0 


4.145 
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E+Dec 


NuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-33.060 


-43.116 


-33.632 


-32.526 


-60.022 


-30.000 


3.0 


3.0 


-12.030 


-17.939 


-12.101 


-11.467 


-19.967 


-7.000 


3.0 


30.0 


1.758 


-7.562 


2.770 


2.162 


1.966 


-0.200 


7.0 


1.0 


-27.240 


-43.116 


-27.812 


-26.706 


-54.202 


-30.000 


7.0 


3.0 


-10.315 


-17.913 


-10.386 


-9.751 


-18.251 


-7.000 


7.0 


30.0 


1.761 


-7.560 


2.774 


2.165 


1.969 


-0.200 


11.0 


1.0 


4.424 


-43.116 


5.531 


5.097 


3.614 


-1.000 


11.0 


3.0 


4.448 


-17.912 


5.556 


5.132 


3.702 


-1.000 


11.0 


30.0 


4.958 


-7.317 


6.169 


5.313 


5.262 


-0.200 
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52 
24 


CR 28 — *■ 


52 MN 27 






Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-29.918 


-30.410 


-29.845 


-30.688 


-67.448 


-40.000 


3.0 


3.0 


-10.505 


-12.380 


-10.326 


-10.872 


-22.257 


-10.000 


3.0 


30.0 


1.363 


-1.420 


2.421 


1.615 


1.753 


-0.300 


7.0 


1.0 


-35.738 


-32.876 


-33.575 


-34.094 


-73.268 


-40.000 


7.0 


3.0 


-12.220 


-12.657 


-11.943 


-12.532 


-23.973 


-10.000 


7.0 


30.0 


1.359 


-1.421 


2.417 


1.611 


1.749 


-0.300 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-50.697 


-48.709 


-48.868 


-49.403 


-62.450 


-10.000 


11.0 


30.0 


-2.499 


-4.653 


-1.442 


-2.246 


-2.113 


-0.300 








52 

24 


CR 28 — * 


ok V 29 






Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-27.351 


-32.503 


-27.895 


-29.087 


-67.464 


-40.000 


3.0 


3.0 


-10.522 


-13.094 


-10.470 


-11.337 


-22.153 


-10.000 


3.0 


30.0 


1.681 


-2.123 


2.731 


1.926 


1.783 


-0.300 


7.0 


1.0 


-21.645 


-32.503 


-22.160 


-23.331 


-61.644 


-40.000 


7.0 


3.0 


-8.858 


-13.086 


-8.807 


-9.659 


-20.438 


-10.000 


7.0 


30.0 


1.685 


-2.122 


2.734 


1.930 


1.787 


-0.300 


11.0 


1.0 


4.486 


-32.503 


5.589 


5.027 


2.896 


-2.000 


11.0 


3.0 


4.486 


-13.086 


5.590 


5.026 


2.899 


-2.000 


11.0 


30.0 


4.784 


-1.980 


6.022 


4.979 


5.024 


-0.200 








52 
25 


MN 27 -> 


52 pp 
9fi rt 26 

ZD J 






Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-17.960 


-22.878 


-18.268 


-33.772 


-59.707 


-40.000 


3.0 


3.0 


-8.324 


-10.148 


-8.149 


-11.668 


-20.004 


-10.000 


3.0 


30.0 


1.074 


-5.435 


2.103 


1.244 


1.627 


-0.100 


7.0 


1.0 


-23.780 


-23.313 


-23.436 


-39.541 


-65.527 


-40.000 


7.0 


3.0 


-10.039 


-10.404 


-9.782 


-13.383 


-21.719 


-10.000 


7.0 


30.0 


1.071 


-5.436 


2.100 


1.241 


1.623 


-0.100 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-48.517 


-48.193 


-48.101 


-51.861 


-60.197 


-10.000 


11.0 


30.0 


-2.790 


-8.897 


-1.762 


-2.618 


-2.239 


-0.100 



MN 27 -> 24 CR 28 



Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-6.289 


-6.218 


-6.058 


-5.418 


-41.039 


-30.000 


3.0 


3.0 


-3.610 


-4.881 


-3.281 


-3.029 


-14.153 


-10.000 


3.0 


30.0 


2.092 


-1.377 


3.178 


2.735 


2.279 


-0.400 


7.0 


1.0 


-2.414 


-6.218 


-2.163 


-1.868 


-35.220 


-30.000 


7.0 


3.0 


-2.124 


-4.865 


-1.760 


-1.564 


-12.439 


-10.000 


7.0 


30.0 


2.096 


-1.376 


3.182 


2.738 


2.283 


-0.400 


11.0 


1.0 


4.630 


-6.218 


5.797 


5.241 


4.399 


-0.300 


11.0 


3.0 


4.644 


-4.865 


5.816 


5.245 


4.436 


-0.300 


11.0 


30.0 


5.050 


-1.248 


6.334 


5.524 


5.458 


-0.200 
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52 
26 



FE 



26 



27 



CO 



Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-77.234 


-100 


-77.183 


-78.493 


-88.784 


-10.000 


3.0 


3.0 


-25.866 


-100 


-25.712 


-26.499 


-29.486 


-3.000 


3.0 


30.0 


-0.080 


-100 


0.912 


-0.787 


0.436 


-0.100 


7.0 


1.0 


-83.054 


-100 


-83.003 


-84.313 


-94.604 


-10.000 


7.0 


3.0 


-27.582 


-100 


-27.427 


-28.215 


-31.201 


-3.000 


7.0 


30.0 


-0.084 


-100 


0.908 


-0.791 


0.433 


-0.100 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-66.059 


-100 


-65.905 


-66.692 


-69.679 


-3.000 


11.0 


30.0 


-3.946 


-100 


-2.954 


-4.653 


-3.430 


-0.100 



52 
26 



FE 



26 



II MN 



27 



Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-6.385 


-4.774 


-5.018 


-5.116 


-33.437 


-30.000 


3.0 


3.0 


-3.955 


-4.780 


-3.558 


-4.114 


-12.095 


-8.000 


3.0 


30.0 


2.088 


-0.361 


3.214 


2.420 


2.330 


-0.400 


7.0 


1.0 


-2.569 


-4.774 


-2.207 


-2.781 


-27.619 


-20.000 


7.0 


3.0 


-2.482 


-4.764 


-2.064 


-2.667 


-10.382 


-7.000 


7.0 


30.0 


2.092 


-0.360 


3.217 


2.423 


2.334 


-0.400 


11.0 


1.0 


4.490 


-4.774 


5.664 


4.920 


4.235 


-0.300 


11.0 


3.0 


4.490 


-4.764 


5.666 


4.921 


4.236 


-0.300 


11.0 


30.0 


4.958 


-0.262 


6.274 


5.146 


5.425 


-0.200 










27 CO 25 


- |NI 


24 





Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-78.338 


-100 


-78.287 


-77.844 


-80.970 


-1.000 


3.0 


3.0 


-26.250 


-100 


-26.099 


-25.956 


-27.607 


-0.400 


3.0 


30.0 


-0.437 


-100 


0.538 


-0.977 


0.141 


-0.040 


7.0 


1.0 


-84.158 


-100 


-84.107 


-83.663 


-86.790 


-1.000 


7.0 


3.0 


-27.966 


-100 


-27.814 


-27.671 


-29.322 


-0.400 


7.0 


30.0 


-0.440 


-100 


0.535 


-0.980 


0.138 


-0.040 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-66.443 


-100 


-66.292 


-66.149 


-67.800 


-0.400 


11.0 


30.0 


-4.303 


-100 


-3.328 


-4.843 


-3.725 


-0.040 








52 
27 


CO 25 - 


52 pp 
26 ht 26 







Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-3.443 


1.040 


1.673 


1.725 


-0.451 


-1.000 


3.0 


3.0 


-1.094 


1.087 


1.723 


1.786 


-0.309 


-1.000 


3.0 


30.0 


2.833 


1.097 


4.041 


3.319 


3.170 


-0.300 


7.0 


1.0 


0.257 


1.040 


1.786 


1.782 


0.113 


-1.000 


7.0 


3.0 


0.324 


1.088 


1.839 


1.842 


0.218 


-1.000 


7.0 


30.0 


2.837 


1.098 


4.045 


3.322 


3.174 


-0.300 


11.0 


1.0 


5.350 


1.040 


6.714 


5.665 


5.783 


-0.200 


11.0 


3.0 


5.386 


1.088 


6.752 


5.724 


5.813 


-0.200 


11.0 


30.0 


5.484 


1.182 


6.894 


5.857 


5.966 


-0.200 
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1 CU 23 




Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


ProEn 


PPEm 


3.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


3.0 


3.0 


-35.406 


-100 


-35.255 


-100 


-35.931 


0.000 


3.0 


30.0 


-0.943 


-100 


0.022 


-100 


-0.484 


0.000 


7.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


7.0 


3.0 


-37.121 


-100 


-36.971 


-100 


-37.646 


0.000 


7.0 


30.0 


-0.946 


-100 


0.019 


-100 


-0.487 


0.000 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 


-100 


11.0 


3.0 


-75.598 


-100 


-75.448 


-100 


-76.124 


0.000 


11.0 


30.0 


-4.809 


-100 


-3.844 


-100 


-4.350 


0.000 








52 
28 


Nl 24 -> 


1? CO 25 






Aden 


T9 


E-Cap 


E+Dec 


NuEn GamEn ProEn PPEm 


3.0 


1.0 


-3.548 


1.264 


1.983 


0.575 - 


0.657 -2.000 


3.0 


3.0 


-1.245 


1.264 


1.985 


0.576 - 


0.493 -2.000 


3.0 


30.0 


3.072 


1.714 


4.276 


1.871 


3.919 -0.090 


7.0 


1.0 


0.153 


1.264 


2.039 


0.602 


0.432 -2.000 


7.0 


3.0 


0.175 


1.264 


2.042 


0.602 


0.506 -2.000 


7.0 


30.0 


3.075 


1.714 


4.279 


1.874 


3.923 -0.090 


11.0 


1.0 


5.363 


1.264 


6.709 


4.004 


6.180 -0.090 


11.0 


3.0 


5.364 


1.264 


6.710 


4.004 


6.181 -0.090 


11.0 


30.0 


5.757 


1.769 


7.151 


4.282 


6.683 -0.050 








53 r A 

20 LA 33 - 


♦ 2ISC32 




Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-5.299 


1.026 


1.636 


1.414 


-2.191 


-3.000 


3.0 


3.0 


-2.272 


1.019 


1.628 


1.429 


-2.130 


-3.000 


3.0 


30.0 


1.847 


0.999 


3.063 


2.255 


2.466 


-0.300 


7.0 


1.0 


-11.119 


1.015 


1.618 


1.403 


-2.512 


-4.000 


7.0 


3.0 


-3.976 


1.009 


1.612 


1.419 


-2.402 


-3.000 


7.0 


30.0 


1.844 


0.998 


3.060 


2.252 


2.462 


-0.300 


11.0 


1.0 


-100 


-74.690 


-75.284 


-74.799 


-100 


-30.000 


11.0 


3.0 


-42.453 


-25.852 


-25.981 


-25.911 


-36.080 


-9.000 


11.0 


30.0 


-2.002 


-1.935 


-0.705 


-1.267 


-1.394 


-0.600 








11 CA 33 - 


53 u 
19 K 34 





Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-89.435 


-100 


-90.017 


-89.537 


-100 


-20.000 


3.0 


3.0 


-30.465 


-100 


-30.562 


-30.528 


-39.408 


-8.000 


3.0 


30.0 


-1.637 


-100 


-0.669 


-1.673 


-2.213 


-1.000 


7.0 


1.0 


-83.615 


-100 


-84.197 


-83.717 


-100 


-20.000 


7.0 


3.0 


-28.750 


-100 


-28.846 


-28.812 


-37.693 


-8.000 


7.0 


30.0 


-1.634 


-100 


-0.665 


-1.669 


-2.209 


-1.000 


11.0 


1.0 


-0.733 


-100 


-0.430 


-1.359 


-1.131 


-0.050 


11.0 


3.0 


-0.236 


-100 


0.388 


-0.336 


-1.058 


-0.400 


11.0 


30.0 


1.952 


-100 


3.028 


1.905 


1.635 


-1.000 
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53 
21 



SC 



32 



53 T| 
22 1 1 



Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-4.709 


1.514 


2.081 


1.810 


-2.224 


-4.000 


3.0 


3.0 


-1.678 


1.495 


2.060 


1.828 


-1.931 


-3.000 


3.0 


30.0 


2.714 


1.495 


3.895 


2.638 


3.435 


-0.100 


7.0 


1.0 


-10.529 


1.501 


2.058 


1.797 


-2.819 


-4.000 


7.0 


3.0 


-3.382 


1.483 


2.039 


1.815 


-2.367 


-4.000 


7.0 


30.0 


2.711 


1.494 


3.892 


2.635 


3.432 


-0.100 


11.0 


1.0 


-100 


-77.703 


-78.296 


-79.700 


-100 


-30.000 


11.0 


3.0 


-41.859 


-27.269 


-27.396 


-27.945 


-36.753 


-9.000 


11.0 


30.0 


-1.136 


-1.550 


0.068 


-0.877 


-0.410 


-0.200 








53 
21 


SC 32 - 


lu CA 33 







Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-55.730 


-100 


-56.307 


-56.594 


-80.171 


-20.000 


3.0 


3.0 


-18.309 


-100 


-18.390 


-18.720 


-27.092 


-8.000 


3.0 


30.0 


1.613 


-100 


2.631 


1.762 


0.161 


-2.000 


7.0 


1.0 


-49.910 


-100 


-50.487 


-50.774 


-74.351 


-20.000 


7.0 


3.0 


-16.594 


-100 


-16.674 


-17.004 


-25.376 


-8.000 


7.0 


30.0 


1.617 


-100 


2.635 


1.765 


0.164 


-2.000 


11.0 


1.0 


1.647 


-100 


2.412 


1.437 


1.985 


-0.200 


11.0 


3.0 


2.183 


-100 


3.126 


2.296 


2.001 


-0.700 


11.0 


30.0 


4.756 


-100 


5.990 


4.900 


3.665 


-1.000 
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23 V 30 







Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-6.382 


-2.012 


-1.645 


-2.087 


-27.002 


-20.000 


3.0 


3.0 


-3.082 


-1.733 


-1.258 


-1.900 


-9.791 


-8.000 


3.0 


30.0 


2.396 


0.449 


3.506 


2.769 


2.736 


-0.200 


7.0 


1.0 


-12.202 


-2.116 


-1.719 


-2.674 


-31.980 


-30.000 


7.0 


3.0 


-4.790 


-1.775 


-1.307 


-2.330 


-11.458 


-9.000 


7.0 


30.0 


2.392 


0.448 


3.502 


2.766 


2.733 


-0.200 


11.0 


1.0 


-100 


-95.742 


-96.339 


-97.574 


-100 


-50.000 


11.0 


3.0 


-43.267 


-32.509 


-32.651 


-33.621 


-49.845 


-20.000 


11.0 


30.0 


-1.460 


-2.627 


-0.349 


-1.061 


-1.125 


-0.200 



53 
22 



Tl 



31 



53 
21 



SC 



32 



Aden 


T9 


E-Cap 


E+Dec 


NuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-49.515 


-100 


-50.092 


-50.844 


-84.119 


-30.000 


3.0 


3.0 


-16.708 


-100 


-16.790 


-17.711 


-27.946 


-10.000 


3.0 


30.0 


1.535 


-100 


2.546 


1.908 


0.863 


-1.000 


7.0 


1.0 


-43.695 


-100 


-44.273 


-45.024 


-78.300 


-30.000 


7.0 


3.0 


-14.992 


-100 


-15.075 


-15.995 


-26.230 


-10.000 


7.0 


30.0 


1.539 


-100 


2.549 


1.911 


0.867 


-1.000 


11.0 


1.0 


3.471 


-100 


4.328 


4.193 


1.740 


-2.000 


11.0 


3.0 


3.473 


-100 


4.334 


4.196 


1.771 


-2.000 


11.0 


30.0 


4.756 


-100 


5.959 


5.099 


4.436 


-0.700 
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53 
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V 
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CR 



Aden 


T9 


E+Cap 


E-Dec 


ANuEn 


GamEn 


NeuEn 


PNEm 


3.0 


1.0 


-6.717 


-1.840 


-1.616 


-2.133 


-33.662 


-30.000 


3.0 


3.0 


-3.629 


-1.888 


-1.649 


-2.123 


-12.149 


-10.000 


3.0 


30.0 


2.037 


0.069 


3.140 


2.313 


2.676 


-0.200 


7.0 


1.0 


-12.537 


-1.953 


-1.773 


-2.326 


-37.635 


-30.000 


7.0 


3.0 


-5.335 


-1.974 


-1.777 


-2.289 


-13.695 


-10.000 


7.0 


30.0 


2.034 


0.068 


3.137 


2.310 


2.673 


-0.200 


11.0 


1.0 


-100 


-100 


-100 


-100 


-100 
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EXPLANATION OF TABLE 



TABLE B. fp-Shell Nuclei Weak Rates of Neutron-Rich and Proton-Rich Nuclei in 
Stellar Matter 



The weak interaction rates for neutron-rich and proton-rich nuclei are calculated using the mass formula of 
Myers and Swiatecki |3(J. The explanation of the table is similar to that of Table A. The superscript (MS) on 
each parent nucleus indicates that the mass formula |3(J is used for the calculation of Q-values and separation 
energies. 55 Zn is an unstable nucleus according to the calculations of |3(J. Hence for the case of 55 Zn the 
calculations of the weak interaction rates are presented using the mass formula of [^9j which is also indicated 
by the superscript (MN) on this nucleus. 
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